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j 要 


孤子 理论 是 非 线性 科学 中 一 个 十 分 重要 的 分 支 ， 它 在 物理 
学 和 其 它 的 许多 领域 中 有 者 越 来 越 广泛 的 应 用 。 而 微 扰 理论 作为 
孤子 理论 中 的 一 个 重要 方面 ， 关 于 它 的 研究 也 日 益 引 起 人 们 的 关 
注 。 到 目前 为 止 ， 人 们 已 经 发 展 了 多 种 有 效 的 微 扰 方法 ， 它 们 大 
体 可 分 为 两 大 类 ;一 类 是 逆 散 射 微 拢 方法 。 它 能 成 功 地 处 理 各 种 
复杂 的 微 扰 问题 ， 在 理论 上 有 着 重要 的 学 术 价值 ， 但 其 思路 较 迁 
回 曲折 ， 数 学 计算 繁琐 。 另 一 种 是 直接 微 扰 方法 。 它 的 最 常用 
的 系统 方法 是 将 孤子 方程 线性 化 后 再 按 Jost 函数 的 平方 作 微 扰 
展开 。 由 于 Jost 解 是 逆 散 射 变换 (IST) 过 程 中 所 出 现 的 一 种 特 
解 ， 所 以 这 种 微 扰 方法 并 没有 完全 脱离 对 IST 的 依赖 ， 与 其 仍 有 
或 多 或 少 的 关联 。 而 且 关 于 直接 微 扰 法 的 含 时 、 空 的 偏 微 分 算 子 
的 本 征 函数 系 ， 前 人 并 没有 找到 直接 求解 它 的 方法 ， 都 是 利用 它 
与 IST 间 的 某 种 联系 来 求 得 的 。 最 近 ， 楼 在 前 人 的 基础 上 发 展 了 
一 种 新 的 直接 微 扰 方 法 ， 它 完全 摆脱 了 对 IST 的 依赖 ， 思 路 和 计 
算 较 为 简便 。 而 且 利用 这 种 直接 微 扰 方 法 得 到 的 微 扰 解 比 用 其 他 
微 扰 方法 得 到 的 解 要 丰富 得 多 。 

本 文 就 楼 的 直接 微 扰 方 法 作 了 较为 全 面 、 系 统 的 研究 ， 主 
要 工作 内 容 如 下 : 1. 从 可 积 的 非 线性 薛 定 刘 方 程 入 手 ， 将 利用 
直接 微 扰 方法 得 到 的 近似 解 与 其 精确 解 进行 了 对 比 、 分 析 ， 验 证 
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了 直接 微 扰 方法 的 可 靠 性 。 微 扰 小 参量 e 越 小 ， 相 应 的 近似 解 与 
精确 解 越 符合 ， 此 时 ， 直 接 微 扰 方法 给 出 的 微 扰 方 程 的 近似 解 也 
就 越 可 靠 。2， 以 不 同 的 微 扰 非 线性 薛 定 刘 方 程 为 实际 例子 ， 例 
dn, SA 3ERRTEBOEISZ EE. (241) VETERE IS f. RA 
述 直接 微 扰 方法 的 具体 运用 和 适用 性 。 通常， 针对 微 扰 非 线性 薛 
定 刘 方 程 的 不 同形 式 的 解 需 要 采用 不 同 的 微 扰 方法 对 其 进行 分 别 
处 理 ， 但 是 ， 从 本 文 的 求解 过 程 中 发 现 ， 采 用 楼 的 直接 微 扰 方 法 
可 以 得 到 多 种 不 同形 式 的 微 扰 解 。 而 且 微 扰 方 程 的 零 阶 解 可 以 为 
无 微 扰 非 线 性 薛 定 刘 方 程 的 任意 一 个 精确 解 ， 它 们 的 一 阶 修正 则 
可 以 是 无 微 扰 方程 的 无 穷 多 个 对 称 的 任意 一 个 。 借 助 符号 处 理 软 
fF Matlab, RIEKA T TREPERE ERES ERER ES 
方程 的 数值 模拟 ， 结 果 表 明 ， 在 e 值 的 较 大 范围 内 ， 近 似 的 解析 
解 与 数值 解 吻合 得 非常 好 。 另 一 方面 ， 非 线性 莅 定 户 方 程 族 的 近 
似 求解 的 实现 也 为 楼 所 发 展 的 直接 微 扰 方 法 提供 了 几 个 重要 的 实 
例 ， 使 直接 微 扰 方法 得 到 进一步 的 完善 和 充实 。3. 基于 直接 微 扰 
方法 给 出 的 可 靠 的 近似 解 ， 我 们 还 给 出 在 微 扰 作用 下 牟 子 的 各 个 
参数 的 解析 表达 式 。 无 论 是 对 于 单 孤 子 还 是 多 匠 子 ， 在 线性 损耗 
的 作用 下 ， 它 们 的 振幅 和 能 量 都 将 随 着 时 间 的 增加 而 指数 衰减 直 
至 零 。 不 过 ， 微 扰 的 影响 并 不 会 改变 抓 子 的 传输 速度 。 特 别 地 ， 
对 于 多 孤子 ， 孤 子 的 碰撞 性 质 也 不 会 受 微 扰 影 响 而 发 生 改变 。 也 
就 是 说 ， 原 来 等 幅 传输 的 弹性 碰撞 的 两 孤子 在 碰撞 后 仍然 保持 以 
相等 的 振幅 传输 ， 它 们 之 间 也 不 发 生 能 量 的 转移 和 交换 。 
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THE DIRECT PERTURBATION METHOD AND 
APPROXIMATE SOLUTIONS OF NONLINEAR 
SCHRODINGER EQUATION HIERARCHY 


ABSTRACT 


Soliton theory is one of the important branches of nonlinear 
science. It has crescent application in many fields. And perturba- 
tion theory is one of the important parts of soliton theory. So the 
study of perturbation theory has been increasingly received exten- 
sive attention of many researchers. So far, many effective pertur- 
bation methods have been established, they can be roughly divided 
into two kinds. One is the inverse scattering perturbation method, 
which has important academic value. But this technique is incon- 
venient to those who are not familiar with the inverse scattering 
transformation(IST). The other is the direct perturbation method 
where the squared Jost solutions are employed as the basis for per- 
turbation expansion after soliton equations have been linearized. 
The method can not completely get rid of the IST because that 
the squared Jost appear as the particular solution in the process 
of the IST. The most important technique in this approach is to 
find eigenfunctions of a linearized operator associated with the lin- 


earized equation. However, the eigenfunctions are actually derived 


Ve 


by making use of some knowledge of IST. Recently, Lou proposed 
a new direct perturbation method. It relies on no tools of IST, 
and it is very concise and easy to understand. The perturbation 
solutions obtained by this perturbation method are more abundant 


than those of other perturbation methods. 


In this thesis, we study systematically Lou's new direct pertur- 
bation method. The contributions of this thesis can be summarized 
as follows: 1. The reliability of the direct perturbation method is 
verified by comparing the approximate solutions of an integrable 
nonlinear Schródinger equation to its exact solutions. For smaller ε 
value, the agreement between the exact and approximate solutions 
is better, which indicates that the approximate solutions obtained 
through the direct perturbation method are more reliable. 2. We 
take different perturbed nonlinear Schródinger equations (i.e., the 
coupled nonlinear Schródinger equation and the (2+1)-dimensional 
nonlinear Schródinger equation) as real examples to show the ap- 
plicability of direct perturbation method. In usual perturbation 
methods, we have to treat different types of solutions of perturbed 
nonlinear Schródinger equations separately, while in our treatment 
all the perturbation solutions can be obtained at the same time. 
Moreover, the zero-order solutions of perturbed equations can be 


any one of the exact solutions of the unperturbed equations, and 
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the first-order corrections can be any one of their symmetries. By 
making use of the computer algebra software Matlab, we simulate 
the coupled nonlinear Schrödinger equations with linear loss nu- 
merically. It turns out that the approximate analytical solutions 
obtained by the direct perturbation method agree well with the nu- 
merical solutions within a wide range of e values. On the other hand, 
our work enriches Lou's direct perturbation method and offers some 
important examples. 3. Based on the reliable approximate solutions 
of perturbed nonlinear Schródinger equations, the analytical form of 
the parameters of solitons can be well expressed. Under the action 
of the perturbation of linear loss, the amplitudes and energy of both 
single-solitons and multi-soliton decay exponentially with time. But 
there is no change in the velocities of solitons. Especially, the per- 
turbations do not change the collision property of multi-soliton, i.e., 
for two soliton with the same amplitude, they will always propagate 
with the same amplitudes both before and after collision and there 


is no exchange and transfer of energy between them. 
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(一 ) 概述 
1 孤立 子 的 由 来 及 发 展 


有 人 认为 非 线 性 科学 是 20 世纪 继 相对 论 、 量 子 力学 之 后 的 “第 三 次 科学 
革命 ”。 一 位 物理 学 家 山 描绘 “相对 论 的 建立 排除 了 对 绝对 空间 和 时 间 的 牛 
顿 幻想 ， 量 子 力 学 的 建立 则 排除 了 对 可 控制 测量 过 程 的 牛顿 迷梦 ， 而 非 线 性 科 
学 的 建立 更 是 排除 了 拉 普 拉 斯 决定 论 的 可 预见 性 的 狂想 ”。 非 线性 科学 正在 成 
为 跨 学 科 的 研究 前 沿 ， 它 主要 针对 各 学 科 中 有 关 非 线性 的 共性 问题 进行 研究 ， 
特别 是 那些 无 法 从 线性 模型 稍 加 修正 就 可 以 解决 的 问题 ， 以 及 它 自 身 理论 发 展 
所 需要 的 概念 和 方法 。 换 言 之 ， 非 线性 科学 是 揭示 各 种 非 线性 现象 的 共性 ， 发 
展 处 理 它们 的 普 适 方法 。 它 的 研究 主体 是 混沌 、 分 形 和 孤立 子 ， 而 且 这 三 者 彼 
此 互相 联系 。 在 本 论文 中 ， 我 们 将 研究 其 中 的 孤立 子 问题 。 

ILF (Soliton) 的 名 称 ，1965 年 由 Zabusky 和 Kruskal!) 提出 ， 但 最 早 
对 孤立 子 这 种 非 线 性 现象 的 研究 可 以 追溯 到 1834 年 ， 英 国 科学 家 罗 梭 (James 
Scott Russell) 对 爱丁堡 到 格拉 斯 哥 的 运河 中 浅水 面 上 水 波 的 观察 。1844 年 ， 
他 在 《英国 科学 促进 协会 第 14 届 会 议 报告 》 这 份 材料 上 发 表 的 “ 论 波动 ”一 
文中 ， 对 观察 到 的 现象 作 了 生动 的 描述 B 

我 正在 观察 着 一 般 船 的 运动 ， 它 被 一 对 马 拉 着 快速 地 沿 着 狭 窗 的 运河 前 进 ， 当 这 
艘 船 突然 停 下 时 ， 被 船 所 带动 的 大 量 河 水 却 并 不 停止 ， 它 们 积聚 在 船 头 周围 剧烈 地 扰 
动 着 。 然 后 突然 离开 船 头 ， 以 很 快 的 速度 向 前 滚动 ， 形 成 一 个 圆 而 光滑 、 轮 廓 分 明 
的 、 巨 大 的 孤立 的 水 堆 ， 沿 着 运河 继续 前 进 ， 在 行进 过 程 中 没有 明显 的 形状 改变 和 速 
度 减 少 。 我 骑 着 马 继续 跟踪 ， 发 现 它 仍 保持 原 有 的 大 约 90 英尺 的 长 度 和 L5 英尺 的 
浪头 ， 并 以 约 每 小 时 8 至 9 英里 的 速度 向 前 滚动 。 然 后 ， 高 度 逐 渐 减 少 ， 约 1. 2X 
里 后 消失 在 河道 的 转弯 处 ……:， 

这 就 是 1834 年 8 月 Russell 第 一 次 观察 到 的 奇特 现象 ， 他 认为 这 种 孤立 的 
波动 是 流体 运动 的 一 个 稳定 解 ， 并 称 它 为 “孤立 波 ”。 斯 科 特 (A.C.Scott)， 
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朱 (F.Y.F.Chu) 和 考 克 劳 林 (D.W.MckLaughlin)) 曾 对 孤立 波 给 出 如 下 的 定 
X: 孤立 波 是 一 种 扰动 ， 它 是 具有 局 域 性 的 行 波 。 孤 立波 到 处 正则 ， 没 有 奇异 
性 ， 而 且 它 不 扩散 ， 与 普通 的 波 不 同 ， 普 通 的 波 包 要 弥散 。10 年 之 后 的 1844 
年 ，Russell 在 浅水 槽 中 做 了 一 些 实 验 ， 用 多 种 方法 激发 ， 看 到 了 相同 的 现象 。 
但 是 Russell 的 学 说 并 未 能 成 功 地 使 当时 的 物理 学 家 信服 ， 对 这 种 波 的 性 质 的 
研究 直到 1895 年 KdV(Korteweg de Vries) 方程 的 提出 和 20 世纪 60 ERR 
射 方法 的 发 展 起 来 后 才 得 到 彻底 解决 。 


关于 孤立 子 ， 目 前 尚 无 公认 的 定义 。 李 政道 向 认为 :在 一 个 场 论 系统 
中 ， 如 果 有 一 个 经 典 解 ， 它 在 任何 时 间 都 束缚 于 空间 的 有 限 区 域内 ， 那 么 ， 这 
样 的 解 就 叫做 经 典 狐 立 子 解 。 可 以 注意 到 ， 这 里 并 没有 给 出 孤子 的 明确 定义 。 
因此 ， 比 较 恰 当 的 做 法 是 把 如 下 三 点 作为 孤子 的 定义 ， 即 ; 

(1) 孤立 子 是 波动 问题 中 的 一 种 能 量 有 限 局 域 解 ; 

(2) 能 在 空间 给 定 区 域 稳定 存在 ; 

(3) 相互 作用 不 改变 各 自 的 特性 。 


李 政 道 等 人 基于 对 基本 粒子 的 研究 ， 又 把 现 有 的 孤立 子 分 为 两 类 ， 拓扑 性 
孤立 子 和 非 拓扑 性 孤立 子 。 前 者 稳定 存在 的 必要 条 件 是 有 简 并 真空 态 ， 即 在 无 
穷 远 处 存在 不 同 的 真空 态 ， 或 者 说 有 不 同 的 边界 条 件 ， 有 孤立 子 解 时 ， 无 穷 远 
处 的 边界 条 件 就 与 没有 孤立 子 解 时 的 不 同 。 非 拓扑 孤立 子 不 需要 简 并 真空 态 ， 
无 论 有 无 孤立 子 ， 在 无 穷 远 处 都 有 相同 的 边界 条 件 。 一 般 地 说 ， 钟 型 分 布 的 
正 、 负 孤 波 及 其 序列 都 是 非 拓扑 的 ， 但 是 kink MAR (其 模 方 或 其 导数 却 是 钟 
型 的 ， 如 光纤 中 基本 了 暗 孤 子 就 是 例子 之 一 ) 是 拓扑 孤子 。 


2 孤立 子 方程 


通常 把 描述 波动 在 空间 内 作 定 域 分 布 ， 而 且 在 有 限时 间 内 表现 相对 稳定 的 
那些 非 线性 方程 ， 都 称 为 孤立 子 方程 。 到 目前 为 止 ， 很 多 物理 和 应 用 数学 领域 
已 经 找到 了 一 大 批 重 要 的 孤立 子 方程 和 方程 组 。 下 面 我 们 列举 一 些 常见 的 具有 
代表 性 的 孤立 子 方程 及 其 对 应 的 孤子 解 : 


(—) ΒΕ -3- 


描写 浅水 波 的 Κάν 方程 ; 
Uy, — Gut, + χες = 0. (1.1) 
KdV 方程 也 是 描述 其 他 许多 领域 中 孤 波 现象 的 模型 ， 如 在 长 波 小 振幅 近似 下 ， 
可 描写 冷 等 离子 体 的 磁 流 体 波 的 运动 ， 非 谱 振 晶 格 的 振动 ;等 离子 体 离子 声 


波 ; 弹性 杆 中 纵向 色散 波 ; 液 、 气 两 种 混合 态 的 压力 波 ， 管 底下 部 流体 底 转 
动 ， 低温 下 非 线性 晶 格 底 声 子 波 包 底 热 激 发 等 。 


1-1: 钟 型 狐 子 
Fig.1-1: Bell soliton 


Tg SPD (大 体形 状 如 图 1-1), Bn 
u(x,t) = -2a^sech?[a(z — 4a?t — το)]. (1:3) 


其 中 a，zo 为 常数 ， 分 别 表征 孤子 的 高 度 和 初始 位 置 。 钟 型 解 表示 波动 定 域 于 
较 小 的 空间 范围 中 ， 一 旦 波动 传 过 后 ， 介 质 即 恢复 原 有 状态 。 

描述 流体 短波 和 小 振幅 极限 的 方程 、 描 述 光 纤 中 孤子 传输 方程 和 描述 光波 
自 聚 焦 的 方程 ， 都 可 妇 结 为 非 线 性 薛 定 刘 (nonlinear Schrodinger) 方程 : 


ity + us, + 2|u[^u = 0. (1.3) 
其 中 , i = VI， 它 有 如 下 形式 的 单 孤 子 解 ; 


u(x,t) = 2asech2a(x — 40t)e-?**9^-9)9. (αγ) (1.4) 
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图 1-2: 包 络 孤子 
Fig.1-2: Envelope soliton 


它 表 示 一 个 被 调制 了 的 平面 波 ， 其 包 络 线 2osech2o(x — 4048) 为 钟 型 珀 子 ， 故 
称 包 络 孤子 ， 形 状 大 体 如 图 1-2. 

描述 相 变 问 题 、 晶 格 缺 陷 、 自 感应 透明 、 自 旋 的 传播 、 铁 磁体 的 Bloch 
墙 、 超 导 Josephson 结 、 超 离子 导体 、 电 荷 密度 、 家 面 外 延 等 方面 的 非 线性 孤 
波 ， 都 可 妇 结 为 正弦 一 戈 登 (sine-Gordon) 方程 


Ust — Ure + sinu = 0. (1.5) 


它 的 解 有 两 种 形式 ， 


图 1-3: (a) 扭 结 解 (b) 反 扭 结 解 
Fig.1-3: (a)kink (b)anti-kink 


(一 ) 概述 «5. 


u(x,t) = 4arctan[et(7vt720/ 1-7]. (1.6) 


AP zo 为 常数 ， 形 状 如 图 1-3。 它 们 表示 在 两 种 具有 相同 能 量 和 不 同 构 形 的 
状态 之 间 的 转变 过 程 ， 当 扭 结 在 空间 某 处 传 过 后 ， 该 处 介质 的 构 形 就 将 发 生 转 
变 。 


图 1-4: 呼吸 子 
Fig.1-4: Breathers 


另 一 种 类 型 是 呼吸 子 解 : 


gsinw(z — t) 


u(x,t) — 4 arctan| ech CEN) 


iF (1.7) 
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RHP +u=1, BRAGA 1-4。 它 表示 一 种 强度 作 周 期 变化 的 波 包 。 
除了 以 上 的 几 种 典型 的 孤子 方程 以 外 ， 人 们 还 推导 出 了 许多 著名 的 具有 孤 
子 解 的 标准 的 非 线性 演化 方程 或 方程 组 ， 如 
修正 Kortewey de Vies(MKdV) 方程 ; 


Ug uus + Use, = 0, (1.8) 
Burgers 方程 ; 
Us F Uus + Uge = 0, (1.9) 
Boussinesq 方程 : 
Ut = Urg + (3u”) oe + Usus, (1.10) 
~ 方程 ; 
Use — Uy = -U H u? (1.11) 


等 等 。 同 时 ， 人 们 也 提出 了 许多 求 标准 的 非 线性 演化 方程 孤子 解 的 方法 ， 如 逆 
散射 方法 (ISTE 7, Hirota 双 线 性 方法 18], Painlevé 分 析 法 910, Darboux 
变换 QI， 对 称 约 化 方法 134, 3E BEA I 3, Bückdund 变换 法 14 等 等 。 
方程 的 解 也 从 孤子 推广 到 反 孤 子 、 呼 吸 子 、 扭 结 解 、 涡 旋 解 、 瞬 子 解 、 磁 单 极 
子 解 等 。 扳 立 子 所 具有 的 广泛 的 实际 应 用 背景 促使 孤立 子 理论 的 研究 和 应 用 均 
得 到 了 迅速 的 发 展 。 


3 微 扰 理 论 的 概述 


标准 的 模型 往往 是 高 度 理想 化 的 ， 在 实际 问题 中 ， 考 虑 到 某 些 实际 因素 ， 
如 阻尼 、 外 加 驱动 等 ， 往 往 要 讨论 包含 修正 项 的 对 应 的 非 线性 偏 微分 方程 ， 这 
时 ， 要 严格 地 求 得 其 精确 解 非常 困难 。 于 是 ， 当 修正 项 〈 即 实际 的 系统 和 这 些 
理想 化 的 模型 之 间 往 往 存 在 的 一 些 细微 差异 ， 我 们 称 之 为 微 扰 ) 可 以 作为 小 量 
时 ， 发 展 实用 的 微 扰 方 法 也 就 成 为 必要 。 由 此 孤子 微 扰 理论 就 应 运 而 生 了 。 由 
于 微 扰 存在 的 普遍 性 ， 研 究 微 扰 对 非 线 性 演化 方程 的 孤子 解 的 影响 的 微 扰 理 
论 也 日 益 成 为 孤子 理论 研究 的 一 个 重要 方面 。 一 般 地 ， 微 扰 可 以 分 为 两 类 : 


C) 概述 .7- 


一 类 是 哈密 顿 微 扰 (Hamitonian Perturbation); 另 一 类 是 耗 散 微 扰 (Dissipative 
Perturbation). 

到 目前 为 止 ， 人 们 已 经 发 展 了 很 多 有 效 的 微 扰 方法 ， 例 如 : () 逆 散 射 微 
扰 方 法 09.021, (i) 哈密 顿 微 扰 方 法 23F-26)， 拉 格 朗 日 微 扰 方法 I, Gic 
种 方法 都 起 源 于 变 分 法 ) (ii) 修正 守恒 律 微 扰 方法 C3-82, Gv) 直接 微 扰 方 
法 93-072, (y) 奇 点 微 扰 方 法 5391-4001, (vi) 线性 微 扰 方 法 4T, 48; (vii)Hirota 
微 扰 方法 9 等 等 。 下 面 着 重 介绍 其 中 的 三 种 微 扰 方 法 。 

(1) 修正 守恒 律 微 扰 方 法 

当 修正 项 很 小 且 远 离 共 振 时 ， 孤 子 的 演化 行为 是 绝热 的 ， 在 这 种 情况 下 ， 
修正 影响 仅 表现 为 孤子 参数 (如 幅度 、 频 率 ) 的 变化 ， 而 不 产生 相应 的 附加 项 。 
绝热 关系 体现 为 非 线性 与 色散 的 平衡 ， 非 线性 效应 与 脉冲 幅度 有 关 ， 而 色散 则 
与 脉冲 宽度 联系 在 一 起 ， 故 绝热 关系 体现 为 (振幅 )x( 脉 宽 )= 常数 。 

下 面 以 微 扰 非 线性 薛 定 廖 方程 为 例 来 说 明 这 种 理论 的 主要 思想 ; 


.0q 1620 - 
izz + τ + lal? q= ieR{q, q Jj? (1.12) 


AP ieRIg,q*] 为 微 扰 项 ，e < 1。 为 描述 孤子 参数 的 变化 ， 可 导出 其 演化 方 
程 。 将 式 (1.12) χη" WER (1.12) MRR xg， 将 所 得 方程 对 T 积分 ， 并 注 
意 到 q,qr ETHLACLA 0, WHERE EET TER 


W= f "aar (1.13) 
和 总 能 量 
W -2Y^n, (1.14) 
i=] 
可 得 
1 (R+ qR')dT (1.15) 
及 - 
zd = je I (GR + arR*)dT (1.16) 
as 


=z f (Gr ο μυς 
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AP W,M, Hy 分 别 为 孤子 的 能 量 、 动 量 与 哈密 顿 量 。 在 微 扰 情 况 下 ， 假 设 方 
程 具有 如 下 形式 的 解 : 


q(T, 2) = 4(Z)sech(n(Z)(T —Ty(Z))e "T Hez, (1.18) 


式 中 (2Z) MoZ 由 下 式 给 出 


(n? — r°). (1.19) 
将 式 (1.18) 代入 (1.15). (1.16), WKF η(Ζ) 5 «(Z) 的 演化 方程 


d e 
3 -5 / (R + qR*)dT, (1.20) 


TEJ GR- ak") 
«让 
27 J οο 
注意 到 ， 式 (1.20) 和 (1.21) 精确 到 eA, RRL, GSM ARKH 
线性 孤子 成 分 ， 但 修正 守恒 律 微 拓 方 法 不 能 对 它们 进行 分 析 。 另 外 ， 象 高 阶 色 
散 等 一 些 属于 哈密 顿 型 的 微 扰 项 ， 它 们 并 不 影响 孤子 幅度 和 中 心 频率 k， 但 
可 能 影响 孤子 的 中 心 位 置 及 相位 。 但 在 修正 守恒 律 方法 中 ， 只 有 孤子 幅度 和 频 
率 两 个 独立 的 变量 ， 难 以 描述 哈密 顿 型 的 微 扰 问题 。 所 以 修正 守恒 律 方法 有 其 
局 限 性 ， 不 能 用 于 分 析 非 孤子 成 分 与 哈密 顿 型 的 微 扰 问题 。 但 在 研究 非 哈 密 顿 
型 微 扰 ， 如 损耗 、 非 线性 增益 、 滤 波 器 等 微 扰 时 ， 修 正 守 恒 律 微 扰 方法 具有 广 
泛 的 应 用 。 在 很 多 情况 下 ， 修 正 守 恒 律 微 扰 方 法 也 适用 于 多 孤子 情况 。 
(2) 逆 散 射 微 扰 方 法 
在 20 世纪 70 年 代 后 期 ，Kaupll5，Karpman 和 Maslovl16! 2-358 1 T jit 
散射 微 扰 方 法 ， 并 在 随后 的 一 些 工作 中 得 以 扩充 17. 18。 关 于 逆 散 射 微 扰 方 
法 的 一 些 早期 的 成 果 ，Karpman 和 Newell 分 别 在 文献 [19] 和 [90] 中 作 了 回 
顾 和 和 总结， 并 且 Abdullaev 和 Khabibullaevi2! 对 此 作 了 补充 。 后 来 ，Kivshar 
和 Malomed?2) 对 这 方面 的 工作 作 了 较为 详尽 的 总 结 。 


(q'R+qR*)dT. (1.21) 


C) ΒΕ -9- 


这 种 方法 的 主要 思想 是 : 运用 Lax 对 形式 ， 孤 子 的 微 扰 方程 形式 为 : 
元 + [L, B] = eR [u]. (1.22) 
此 时 ， 算 子 L 依然 满足 Sturm-Liouvill 方程 
LW = do, (1.23) 


但 是 ， 和 未 受 微 扰 影 响 的 情况 不 一 样 ， 和 已 不 再 独立 于 时 间 了 ， 本 征 函数 ὅ 
也 不 满足 演化 方程 D = BS. ΠΗ͂ Jost 解 的 引入 ， 它 们 的 解析 性 和 渐进 行为 
都 与 不 含 修正 项 的 非 线性 演化 方程 一 样 ， 逆 散射 方程 也 保留 原来 同样 的 形式 ， 
只 是 散射 数据 S(t) = {r(t, 入 ,cn(t), An} 随时 间 的 演化 需 另行 设法 求 出 。 

方程 (1.22) 对 时 间 t 微分 : 


对 于 分 立 谱 得 : 

(E ~ A(T: ~ BV) = Xv -- en[ujv. (1.24) 
对 于 连续 谱 得 : 

全 一 ye - BY,) = —eR[u]V,. (1.25) 


(1.24) 和 (1.25) 是 非 齐 次 常 微分 方程 ， 要 先 求解 其 对 应 的 齐 次 方程 ， 然 后 
利用 参数 变易 法 求解 ， 可 导出 离散 谱 ， 反 射 系数 ， 正 归 化 系数 随时 间 的 演化 关 
系 。 由 (1.24) 和 (1.25) 可 定 出 得 散射 数据 随时 间 的 演化 关系 代入 逆 散 射 方程 定 
出 解 ， 此 解 将 满足 含 修正 项 的 非 线性 演化 方程 。 所 以 整个 方法 的 关键 在 于 : E 
样 确定 散射 数据 随时 间 的 演化 。 

此 方法 的 优点 与 不 足 ，(i) 处 理 孤 子 微 扰 问 题 的 能 力 很 强 ， 能 成 功 地 处 理 很 
多 的 孤子 微 扰 问题 ， 如 KdV, SG 和 NLS 等 等 非 线性 演化 方程 的 孤子 微 扰 问 
题 。 

(i) 思路 迁 回 曲 折 ， 而 且 该 方法 建立 在 IST 的 基础 上 ， 它 要 求 未 受 微 扰 影 
响 的 非 线性 演化 方程 能 够 用 IST 精确 求解 。 因 而 只 适用 于 可 积 系 统 。 而 就 一 般 
情况 而 言 ， 拉 克 斯 (Lax) 对 算 子 较 难 求 得 ， 这 就 使 得 此 方法 的 应 用 有 很 大 的 局 
限 性 。 且 对 于 那些 不 懂得 IST 法 的 人 来 说 ， 想 运用 此 方法 非常 困难 。 
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(3) 直接 微 扰 方法 

此 方法 由 Ostrovskii 和 他 的 同事 创立 B3, 3。 随后 ， 此 方法 在 不 同 的 微 
扰 模 型 的 应 用 中 得 以 发 展 和 完善 : 如 微 扰 的 sineGordon FE 95-38, do 
的 nonliear Schrödinger 方程 89， 微 扰 的 nonlinear Klein-Gordon 方程 40 等 
等 。 

该 方法 的 主要 思想 是 ， 将 微 扰 方法 的 解 在 未 受 微 扰 影 响 的 方程 的 孤子 解 附 
近 展 开 而 导出 一 个 关于 修正 解 的 非 齐 次 线性 化 偏 微分 方程 ， 然 后 利用 与 之 相应 
的 线性 偏 微分 算 子 的 本 征 函数 构造 一 套 正 交 完 备 化 基 矢 来 求解 上 述 偏 微分 方 
程 。 在 直接 孤子 微 拢 方法 的 发 展 中 ，Herman 做 了 很 杰出 的 成 绩 。 他 的 方法 简 
述 如 下 ; 

A. 微 扰 方 程 的 线性 化 

引进 多 重 尺度 ， 将 解 作 渐 进 级 数 展开 u= wo + eu +, 得 到 一 修正 方程 


In, = Fy. (1.26) 
其 中 了 是 一 个 含 时 、 空 变量 的 偏 微分 算 子 。 


B. 正 交 完备 化 基 矢 的 构造 
a) 构造 一 正 交 完备 化 的 基 矢 去 求解 方程 (1.26) 


iQ = dN, (1.27) 
LAQA = vA, (1.28) 


其 中 I^ 是 2 的 伴随 算 子 。 由 逆 散 射 理论 中 的 Jost 函数 的 平方 函数 与 算 子 上 
和 L^ 本 征 函 数 的 关系 来 求解 和 14 的 本 征 函数 系 ， 且 使 其 完备 化 。 
b) 由 线性 算 子 的 反 演 求 得 微 扰 解 


Lu =A, (1.29) 
es / f dA, (1.30) 
iw ji f dA. (1.31) 


C. 微 扰 对 孤子 的 影响 的 分 析 


(一 ) 概述 11- 


消除 ul 中 的 久 期 项 获得 孤子 参数 随时 间 的 演化 ， 一 般 相应 于 线性 算 子 的 
RAS. 

此 方法 的 优 缺 点 。，(i) 此 方法 只 需 知道 无 微 扰 方程 的 精确 的 孤子 解 和 已 线性 
化 方程 算 子 L 的 本 征 函 数 ， 它 不 要 求 无 微 扰 方 程 是 精确 可 积 的 ， 因 此 它 既 适用 
于 可 积 系 统 ， 又 适用 于 非 可 积 系 统 。 

(i) 此 方法 思想 表面 上 看 来 是 与 IST 毫 不 相干 ， 可 实际 操作 过 程 中 巧妙 地 
运用 了 一 些 IST 的 结论 。 所 以 此 微 扰 法 并 没有 完全 脱离 对 IST 的 依赖 ， 与 其 
仍 有 或 多 或 少 的 关联 。 且 关于 直接 法 的 含 时 、 空 的 偏 微分 算 子 的 本 征 函 数 系 ， 
前 人 并 没有 找到 直接 求解 它 的 方法 ， 都 是 利用 它 与 IST 间 的 某 种 联系 来 求 得 ， 
因此 它 间接 地 依赖 于 IST. 

(ii) 直接 微 扰 方 法 只 适用 于 单 孤子 微 扰 问 题 。 

由 于 以 上 的 直接 微 扰 方 法 对 IST 的 间接 依赖 ， 很 大 程度 上 限制 了 它 的 发 
E. HB) 20 世纪 90 ER, MRE 542 发 展 了 一 种 新 的 直接 微 扰 理 论 。 这 
种 方法 完全 摆脱 了 对 IST 的 依赖 ， 具 有 较 大 的 普 适 性 和 实用 性 。 

此 微 扰 方法 由 三 部 分 组 成 : 

(i) 运用 经 典 微 扰 理 论 中 的 微分 展开 法 线性 化 微 扰 的 孤子 方程 ; 

(ii) 求 被 线性 化 的 方程 中 的 线性 算 子 及 其 伴随 算 子 的 本 征 值 问 题 ， 构 造 两 
个 正 交 完备 化 的 本 征 函 数 集 ; 

(ui) 微 扰 对 孤子 的 影响 的 分 析 :， 以 本 征 函 数 集 为 基 来 求解 被 线性 化 的 方 
程 ， 求 得 一 阶 近似 下 的 孤子 解 ， 通 过 消去 久 期 项 的 条 件 求 得 孤子 参数 随时 间 和 
空间 慢 变量 的 变化 情况 。 


在 前 人 的 基础 上 ， 楼 50) 发 展 了 一 种 新 的 直接 微 扰 方法 。 它 完全 摆脱 了 
对 IST 的 依赖 ， 且 思路 直接 ， 容 易 理解 和 接受 ， 实 际 操作 过 程 比 其 他 的 微 扰 方 
法 也 要 简单 得 多 。 更 值得 一 提 的 是 ， 在 利用 此 法 来 处 理 微 扰 问 题 的 过 程 中 ， 我 
们 无 需求 任何 函数 的 本 征 值 ， 并 且 能 很 容易 得 到 微 扰 方 程 的 一 阶 修正 。 而 且 利 
用 这 种 直接 微 扰 方 法 得 到 的 微 扰 解 远 比 用 其 他 微 扰 方 法 得 到 的 解 要 丰富 得 多 。 
此 方法 主要 由 以 下 三 个 部 分 组 成 ; 
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(i) 将 拟 设 的 解 的 形式 代入 微 扰 方 程 ， 得 到 关于 e 的 各 级 近似 方程 ， 

(i) 令 微 扰 小 参量 。 HHAKENRRSTE: 其 中 使 得 零 阶 解 为 未 受 微 
扰 影 响 的 非 线性 演化 方程 的 精确 解 ， 一 阶 修正 为 该 演化 方程 的 对 称 。 

(iii) 解 满足 以 上 情况 时 的 相应 的 相 容 性 条 件 。 

目前 这 种 方法 尚 在 进一步 发 展 完善 之 中 。 


(=) 论文 的 研究 意义 及 其 框架 


非 线性 薛 定 刘 方 程 是 现代 科学 中 最 重要 也 是 最 普遍 的 非 线性 模型 之 一 。 作 
为 描述 波 包 在 弱 非 线性 色散 介质 中 传播 的 普遍 方程 ， 它 出 现在 物理 和 应 用 数学 
的 各 个 领域 中 。 例 如 ， 光 纤 孤 立 子 通信 ， 非 线性 晶 格 的 振动 问题 ， 水 表面 波 问 
题 ， 充 满 流体 的 弹性 管道 中 的 波动 问题 ， 血 管 中 血 流 非 线性 波 问 题 等 都 由 非 线 
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讨论 包含 修正 项 的 方程 ， 这 时 ， 要 严格 地 求解 一 般 是 不 可 能 的 。 于 是 ， 当 修正 
项 可 以 当 作 微 扰 时 ， 发 展 实用 的 微 扰 方 法 成 为 必要 。 由 此 微 扰 理 论 就 应 运 而 生 
了 。 而 且 由 于 微 扰 存 在 的 普遍 性 ， 研 究 微 扰 对 非 线性 演化 方程 的 孤子 解 的 影响 
具有 更 为 实际 的 意义 。 

直到 目前 为 止 ， 为 解决 非 线 性 薛 定 协 方程 的 微 扰 问题 ， 人 们 已 经 提出 了 很 
多 有 效 的 微 扰 方法 ， 包 括 修正 守恒 律 微 扰 方法 ， 逆 散射 微 扰 方法 ， 直 接 微 扰 
方法 ; 奇 点 微 扰 方 法 等 等 。 我 们 知道 ， 其 中 最 著名 的 是 逆 散 射 微 扰 方 法 。 它 
处 理 微 扰 问题 的 能 力 很 强 ， 能 成 功 地 处 理 很 多 复杂 的 孤子 微 扰 问 题 ， 但 其 思 
BIL ST. WF RLS IST 法 的 人 来 说 想 运 用 此 方法 是 非常 困难 的 ， 
另外 ，Ostrovsky 等 人 提出 的 直接 微 扰 方 法 在 实际 操作 过 程 中 巧妙 地 运用 了 一 
HIST 的 结论 ， 故 它 并 未 摆脱 对 IST 的 依赖 。 最 近 ， 在 前 人 的 基础 上 ， 楼 发 
展 了 一 种 新 的 直接 微 扰 方法 。 它 完全 摆脱 了 对 IST 的 依赖 ， 且 思路 简洁 ， 操 作 
简单 ， 具 有 较 大 的 普 适 性 与 实用 性 。 本 学 位 论文 利用 直接 微 扰 方 法 来 实现 对 微 
扰 非 线性 薛 定 廖 方程 的 近似 求解 ， 为 进一步 讨论 方程 的 解 的 性 质 提供 了 很 好 的 
基础 ， 也 为 进一步 讨论 方程 的 形态 提供 了 一 种 很 好 的 工具 。 
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全 文 共 分 为 四 章 : 第 一 章 简要 地 介绍 了 孤立 子 的 由 来 以 及 它 的 发 展 历程 ; 
并 列举 了 几 种 具有 代表 性 意义 的 孤立 子 方程 以 及 它们 的 孤子 解 ， 同 时 也 简要 地 
回顾 了 微 扰 理 论 的 形成 和 发 展 。 第 二 章 从 标准 的 非 线性 薛 定 翰 方程 出 发 ， 介 绍 
了 它 的 研究 情况 及 丰富 的 孤立 子 解 。 简 要 回顾 了 非 线性 薛 定 刘 方 程 的 微 扰 理 论 
的 发 展 历程 并 详细 介绍 了 本 论文 着 重要 运用 的 直接 微 扰 方法 ， 进 而 通过 将 这 种 
直接 微 扰 方法 运用 到 一 个 可 积 的 含 修正 项 的 非 线性 薛 定 廖 方程 来 验证 直接 微 扰 
方法 的 可 车 性 。 论 文 的 第 三 章 ， 将 楼 的 直接 微 扰 方法 应 用 到 各 类 不 同 的 微 扰 非 
ETERS AE, QM NAS IRER ESTEM (2+1) HARE ES 
方程 。 基 于 微 扰 方法 给 出 的 近似 解 ， 我 们 还 分 析 了 微 扰 项 对 孤子 的 各 个 参数 的 
影响 ， 如 孤子 的 振幅 、 传 播 速度 、 相 位 和 波 数 等 等 。 进 一 步 地 ， 利 用 数值 方法 
一 分 步 傅立叶 方法 对 售 线 性 损耗 项 的 看 合 非 线 性 醉 定 记 方 程 进行 了 数值 模拟 。 
通过 数值 解 与 近似 的 解析 解 的 比较 ， 很 好 地 说 明了 直接 微 扰 方 法 的 准确 性 和 有 
效 性 。 最 后 是 总 结 和 展望 。 
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二 直接 微 扰 方 法 


(一 ) FRR EAR 
AE HEEE EIS TE (NLSE) 


iA; + Azz + 2ΙΑΙΑ πα 0, (2.1) 


是 现代 科学 中 最 重要 也 是 最 普遍 的 非 线 性 模型 之 一 。 方程 (2.1) 中 ，4 是 变 
量 上 与 z 的 复 函 数 ，|A4j? = 44*。 作 为 描述 波 包 在 弱 非 线性 色散 介质 中 传播 
的 普遍 方程 ，NLSE 出 现在 物理 和 其 他 的 许多 分 支 中 ， 包 括 非 线性 重子 场 理 
论 5D-631， 凝 聚 态 物 理 5 和， 等 离子 体 物理 55-657， 非 线性 光学 1981-161] 和 
量子 电子 学 并， 流体 物理 (63, 6 和 ， 微 扰 和 相 变 理论 65. 666， 生 物 物 理 (67) 
等 。 正 是 因为 NLSE 所 具有 的 广泛 的 应 用 性 和 重要 性 ， 一 直 以 来 都 吸引 着 一 大 
批 数 学 家 和 物理 学 家 们 的 注意 力 。NLSE 的 最 著名 的 解 是 孤子 解 。 孤 子 的 典型 
特性 时 以 一 个 局 域 波 的 形式 在 与 另 一 个 孤子 相互 作用 的 时 候 ， 保 持 了 类 粒子 的 
性 质 。 特 别 地 ， 在 非 线 性 光学 方面 ， 由 于 近年 来 在 无 线 电 通讯 领域 ， 信 息 高 速 
公路 已 成 为 一 个 热门 话题 。 孤 子 通讯 系统 由 于 其 具有 符合 高 速 通讯 系统 所 要 求 
的 特性 而 将 成 为 信息 高 速 公路 所 需 传 输 媒 介 的 主要 候选 。 光 孤子 由 于 其 特有 的 
性 质 成 为 了 长 距离 、 高 比特 率 光 通 讯 系统 的 信息 载体 。 由 于 在 通讯 和 超 快 信号 
发 送 系 统 中 的 广泛 应 用 ， 超 短 光 脉 冲 在 光纤 中 的 传输 情况 受到 了 极 大 的 关注 。 
最 近 对 全 光 孤 子 通 讯 系 统 的 研究 表明 ， 孤 子 通 讯 可 以 克服 线性 通讯 系统 速度 、 
距离 的 限制 。 众 所 周知 ， 皮 秒 光 孤子 的 传输 正 是 由 标准 的 NLSE 来 控制 。 所 
以 ， 在 预言 孤子 存在 的 可 能 性 ， 以 及 Mollenauer, Stolen 和 Gordon 的 实验 发 
现 之 后 ，NLSE 的 光 孤 子 被 认为 是 自然 数据 比特 和 超 高 速 光 学 通讯 系统 的 新 生 
代 重 要 选择 。 

在 Gardner, Greene, Kruskal 和 Miural68| 于 1967 年 为 求解 KdV 方程 
创立 了 IST 之 后 ，Zakharov 和 Shabad® 于 1971 年 将 它 推广 到 求解 NLSE 
FF. Zakharov 和 Shabad 建立 的 IST 是 求解 NLSE 和 其 他 类 似 的 方程 的 最 系 
统 的 方法 。 并 在 无 反射 的 情况 下 ， 导 出 了 决定 NLSE 多 孤子 解 的 线 代数 方程 


C) AFR EAB -15. 


组 1691, 
NLSE 的 最 简单 的 单 孤 子 解 可 以 写成 : 
A(z, t) = —216(z, t)sech(2I(z — zo + 4kt)}. (2.2) 
其 中 ， 
G(x, t) = ee tt Pio, (2.3) 


这 里 ，k, 1, do 和 αρ 分 别 为 表征 孤子 速度 、 振 幅 、 初 始 相位 和 初始 位 置 的 四 个 
实 常 数 。 


RIM FAFA: " 
us 
Ap 22 D (24) 
其 中 ， 
Dz = |(m -- με) 3 νὰ + νᾶ] 
424% [tanh 0, tanh 0, 一 gechbgisechb。 cos(q1 — φο)], (2.5) 


No = [μι -- p2)? + vd — νὰ] — ἵΏναίμι — μι) tanh 62}1,8ech0,e- 

+{[(p — He)? — v? + v2] + 1214 (ux — μο) tanh 6, }vasech6.e—”, (2.6) 

随后 ，Hirota 用 以 双 线 性 微分 方程 为 基础 的 直接 法 ， 给 出 了 NLSE 的 孤子 
解 的 显 式 。 假 设 方程 的 (2.1) 的 解 可 家 示 为 : 

υ-ὖ (2.7) 

并 且 不 失 一 般 性 ， 分 母 Sf 设 为 实 函 数 。 另 设 实 函数 f TET MAKER 
级 数 ， 而 复 函 数 9 KT 的 奇 次 寒 展 开 级 数 ， 

f(z.t) =14 fg 4 fei x... y pong uuu, (2.8) 

g(z,t) = ge 4- GPA +. p g@i Den 4 .... (2.9) 


将 (2.7)-(2.9) 代入 方程 (2.1)， 并 令 T 的 同 次 宕 系数 相等 得 NLSE 的 单 孤子 
jg (70), 


ο. 
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n=e™. 
双 孤 子 解 为 ; 


2 = en 十 ez 十 gi moni 012819623 + enormen) +6124 8144804 ， 


fo=1t entnitéis + etn this 十 em Γη +8343 θαι 
peni +13 812-813-0144 025-024 +034 i 
其 中 
eeu = 2(k 一 Κρ)’, 
ef = NR -- hy), 
1 
Ma (j,1 = 1,2 
nj = WT + kt + nf, Wi = ik. (j = 1,2) 
这 里 k; (j=1,2) 均 为 任意 复 常数 。 


用 Hirota 方法 得 到 的 τι 孤子 解 所 对 应 的 分 母 与 分 子 可 表示 成 ; 


2n 2n 
ΓΞ D Ai (peri μη DIL mb 


p=0,1 
gn = Ao(peXsar Mmi Y Hy pry 
p=0,1 
其 中 
Mmt = (7 =1,2,...,n) 
eie = ΠΕΤΕΡ (j,1=1,2,..4n) 


efit = 2(k; 一 ky, (1 <l= 2,3, suey) 


εσυ = Uk k =e, (j <1 = 2,3, ann) 


(2.10) 


(2.11) 


(2.12) 


(2.13) 
(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


(2.19) 
(2.20) 
(2.21) 


(2.22) 


Ai(u), Axu) 表示 当 ml =1,2,...,.n) 取 所 有 可 能 的 0 或 1 时 ， 还 需 分 别 满足 


条 件 : 


Dy = usns 


j=1 251 


(2.23) 
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STONES (2.24) 
j=1 1-1 
另外 还 有 其 他 的 各 种 求解 非 线性 偏 微分 方程 的 方法 也 被 用 来 求解 NLSE, 
如 Darboux 变换 ，Backlund 变换 ，Painlevé 分 析 法 ， 截 断 展开 法 等 等 。 这 里 
就 不 再 一 一 详细 介绍 了 。 

在 实际 应 用 中 ， 诸 如 非 线 性 醉 定 坦 孤子 的 牛顿 粒子 行为 、 静 孤子 现象 、 导 
周期 外 力 驱 动 下 出 现 的 分 形 以 及 光纤 的 损耗 等 各 种 行为 对 孤子 的 影响 引起 了 研 
究 者 的 广泛 兴趣 。 而 要 考虑 这 种 种 的 影响 ， 就 需要 在 基本 方程 上 添加 不 同 的 修 
正 项 。 这 时 ， 我 们 就 需要 处 理 包 含 修正 项 的 NLSE。 对 于 这 样 的 方程 ， 如 果 修 
正 项 可 以 当 作 小 量 ， 即 能 看 作 微 扰 ， 那 么 发 展 有 效 的 微 扰 方法 就 成 为 必要 ， 而 
这 些微 扰 方 法 也 正 是 处 理 实际 问题 所 不 可 缺少 的 。 

Kaup!!5! 在 研究 NLSE 受阻 尼 影 响 时 ， 通 过 将 包含 微 拓 的 动力 学 方程 转换 
到 散射 空间 ， 首 先 建立 了 NLSE 的 微 扰 理 论 。 随 后 ，Keener 和 Mclaughlin!!! 
使 用 格林 函数 法 ， 以 及 Ichikawal39 研究 Langmuir 波 包 受 非 线性 Landau 阻尼 
影响 时 采用 的 “绝热 近似 法 ”都 对 该 方程 的 微 扰 理论 作出 了 重要 的 贡献 。 其 
中 Keener 和 Mclaughlin 利用 双重 尺度 技巧 和 格林 函数 建立 了 一 套 形式 上 独立 
于 IST 的 直接 微 扰 方 法 。 但 我 们 明显 地 可 以 看 出 ， 展 开 基 矢 是 依靠 IST RB 
的 。Kaup 后 来 的 工作 、 以 及 Herman. 73) 发 展 的 直接 微 扰 方法 也 存在 着 同 
样 的 问题 ， 它 们 实际 是 类 似 的 。Herman 使 用 的 展开 基 矢 是 利用 算 子 的 本 征 函 
数 和 IST 里 面 的 Jost 函数 的 关系 ， 以 及 Kaup 的 工作 作为 基础 来 构建 的 。 这 
些 基 矢 的 构造 包含 了 相当 复杂 的 推导 。 之 后 ， 颜 4L 42 提出 的 直接 微 扰 方 法 
虽然 改善 了 对 IST 的 依赖 关系 ， 操 作 起 来 也 更 简单 。 但 是 ， 就 多 孤子 微 扰 问 
题 ， 这 种 直接 方法 也 还 是 没 能 很 好 地 给 以 解决 。 所 以 在 很 大 程度 上 ， 这 种 直接 
微 扰 方 法 还 是 限制 了 人 们 对 某 些 实际 问题 的 处 理 。 

在 随后 的 章节 中 ， 我 们 将 采用 楼 提出 的 直接 微 扰 方 法 来 研究 各 类 NLSE 的 
微 扰 问题 ， 楼 的 方法 不 仅 摆脱 了 对 逆 散 射 方法 的 依赖 ， 使 得 实际 的 操作 过 程 大 
大 简化 ， 更 重要 的 是 这 种 直接 微 扰 方法 能 够 同时 得 到 各 种 不 同类 型 的 微 扰 解 ， 
包括 单 孤 子 的 、 多 孤子 的 以 及 周期 性 的 解 等 等 。 


(18. 浙江 师范 大 学 硕士 论文 : 二 直接 微 扰 方 法 


(二 ) 直接 微 扰 方法 
这 一 节 中 ， 我 们 将 对 楼 在 文献 [50] 中 提出 的 直接 微 扰 方 法 作 详细 的 氢 述 。 
并 就 该 方法 的 可 靠 性 进行 简单 分 析 。 


光纤 损耗 是 描述 光 信 号 在 光纤 中 传输 时 能 量 损失 的 一 个 参数 。 光 纤 每 单位 
长 度 的 损耗 ， 直 接 关系 到 光纤 通信 系统 传输 距离 的 长 短 。 形 成 光纤 损耗 的 原因 
很 多 ， 有 来 自 光纤 本 身 的 损耗 ， 也 有 光纤 与 光源 的 看 合 损耗 以 及 光纤 之 间 的 连 
接 损耗 。 考 虑 损耗 后 的 NLSE 具有 如 下 形式 [74; 


iA, + Azz 1-2|Α|3Α = -ic4. (2.25) 
RP, e= ?78/lio|， 
当 满 足 条 件 
c£ -4:«1 (2.26) 


Hf. 31 (225) 右边 的 项 可 当 作 微 扰 项 ， 这 就 是 我 们 接 下 来 所 要 研究 的 微 
扰 NLSE。 

由 于 考虑 到 微 扰 方 程 的 零 阶 解 应 该 为 无 微 扰 NLSE 的 精确 解 ， 我 们 将 微 
扰 NLSE(2.25) 中 的 场 量 4 展开 成 : 


A = ef AE, τ, ε) = ere DEA C τὴ + eA (£, τ) + O(e)]. (2.27) 


其 中 
a = a(z,t), be b(z,t), €=€(z,t,€), τ Ξτία, ἐ, ε) (2.28) 


均 为 给 定 变 量 的 实 函 数 。 这 里 ，&, r 需 满足 以 下 性 质 : 
{é,7} 3 {x,t}. (2.29) 


也 就 是 说 ， 当 e 一 0 Rf, Ao 为 标准 NLSE 的 精确 解 。 将 展开 式 (2.27) RAX 
程 (2.25)， 得 到 如 下 近似 方程 ; 


Aoeeéz + 26" | Ao] Ao + iAorTs + ε(ἑΑ1ττι + Argel? + 4€%|Ag|2Ar 


+26 AGAT) + ε(ίαι — be) Ao + €(Aze + ibas) Ao + ieAo + Exe Ang 


(=) 直接 微 扰 方 法 -19- 


"Fi& Ao + 2εξε(ας + tbe) Age + 26,7, Age, + Aor Tz, = Ο(ε)). (2.30) 
由 于 考虑 到 Er 有 如 (2.29) 式 的 性 质 ， 所 以 Er te, bee 1928 e 的 一 阶 项 。 今 6 
的 相同 次 器 的 系数 等 于 零 ， 得 到 如 下 各 级 近似 方程 : 
Aog? + 2εῦέα| Aol?2 A + iAoni = 0, (2.31) 
ε(ἰΑιττι + Aie + 4e? | Ag]? Ay T 26** 4245) 
Γείίαι ~ ὃν) Ao 十 e(azz + iba.) Ao + icAo + Erz Aog 
Fit; Age + 2€Ex (az + ib.) Aog t 2£»Tz Ager + Αρττω = 0. (2.32) 
由 于 方程 (2.31) 中 的 场 量 A 是 不 显 含 6 的 ， 所 以 ET 与 a 的 依赖 关系 可 以 表 
示 成 ; 
& =e", n=. (2.33) 
XXFÉ, Αρ 即 为 以 下 NLSE 
i4or 十 4ote 十 24o24=0 (2.34) 


的 精确 解 。 注 意 到 ， 此 时 NLSE(2.34) 以 &,7 为 变量 。 另 外 ， 从 展开 式 (2.27) 
中 不 难 发 现 ， 一 阶 修正 Αι 正 是 NLSE(2.31) 的 对 称 ， 也 就 是 说 ，4; BF 
程 (2.34) 的 线性 化 方程 


Aye + Arge + 4| Ao Αι + 242A] = 0 (2.35) 


的 解 。 我 们 知道 ，NLSE 的 对 称 已 经 被 很 多 的 物理 学 家 和 数学 家 们 研究 
过 /8。 它 的 最 简单 的 两 个 形式 可 以 写成 ， 


Αι = Aoc, A; = Aor- (2.36) 


对 于 给 定 的 一 个 非 平庸 解 A 1 (2.32), (2.35) 和 (2.33) 在 任意 的 空间 
© 处 和 时 间 ¢ 下 都 是 相 容 的 ， 这 就 要 求 方程 (2.32) 中 其 余 的 项 均 为 等 ， 即 


“(i+ Eaa) "2, + ibe =O, 


ilar + ibi) + (Gee + idee) +i=0, τ, -- 0. (2.37) 
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通过 简单 的 计算 得 到 
a = —2t, b = T 


1 — 74t 
= get r= n, 2.38 
E= re”, τ dc (2.38) 


于 是 ， 将 函数 a, b 和 变量 &,7 代 回 展开 式 (2.27) 中 ， 就 可 以 确定 微 扰 NLSE(2.25) 
的 近似 解 。 

通过 以 上 对 直接 微 扰 方 法 的 介绍 ， 我 们 不 难 发 现 ， 微 扰 解 中 的 ho 可 以 
是 NLSE(2.31) 的 任意 一 个 解 ， 包 括 单 孤子 解 ， 多 孤子 解 ， 甚 至 呼吸 子 解 等 。 
它们 的 具体 形式 可 以 参见 本 章 第 一 节 。 其 中 Αι 则 可 以 为 方程 (2.31) HEH 
多 个 对 称 中 的 任意 一 个 。 例 如 ， 同 样 的 零 阶 解 4o 可 以 跟 不 同 的 一 阶 修正 Αι 
组 成 不 同 的 近似 解 ， 而 且 ， 不 同 的 零 阶 解 Αρ 也 可 以 有 相同 的 Αι 作为 一 阶 修 
正 。 所 以 ， 采 用 楼 的 直接 微 扰 方 法 得 到 的 微 扰 解 比 其 他 微 扰 方法 得 到 的 解 要 让 
富 得 多 。 

ERAH NLSE(2.25) 的 零 级 解 ho 为 单 孤 子 解 (2.2)-(2.3)， 则 相应 的 微 扰 
解 可 表示 成 ; 


—2lk + 2lke*t 


A = —2le~*esech| — 2lze^?* + οἱαρ] 


E —4et 
x{1 — 2ick + 2el tanpe +e — 2lxe~** + 2]σρ]}. (2.39) 
其 中 
A = Aot; 
2 (a1. ρ-ε 
y= + eee — 2kze"?* — ġo. (2.40) 


前 面 我 们 提 到 ， 人 们 已 经 利用 各 种 微 扰 方 法 对 方程 (2.25) 作 过 研究 ， 也 得 
到 了 各 种 不 同形 式 的 微 扰 解 。 例 如 (以 单 孤子 解 为 例 ) 
(1) 逆 散射 微 扰 方法 
将 微 扰 解 设 成 ; 
A = 2lsechZe-y， (2.41) 


(=) 直接 微 扰 方法 -21- 


Z = 2l(z -- €), (2.42) 
q(z,t) = 2k(z — 6) + ô. (2.43) 
通过 计算 后 得 到 其 中 各 参量 如 下 À, 


k = ko, (2.44) 

€ = £o — 4kot, (2.45) 
l2le?*5, , (2.46) 

5 = by — Akt — Ba — e7), (2.47) 


以 上 的 量 凡 是 带 下 标 0 的 都 是 任意 常数 (FRI). 34 6— 0, 
5 = by — 4(k + 2t, (2.48) 


l2 d. (2.49) 


在 这 样 的 极限 下 ， 微 扰 解 退化 为 无 微 扰 的 NLSE 的 单 孤子 解 。 

(2) 颜 的 直接 微 扰 方法 

利用 该 方法 得 到 的 k, &, 1 的 形式 与 (2.44)-(2.46) 完全 相同 ， 唯 一 不 同 的 
是 [4 

δ = by — (4ἱῷ — e)t + See. (2.50) 

但 是 ， 当 e 一 0， 微 扰 解 却 不 能 退化 为 无 微 扰 的 NLSE 的 单 孤子 解 。 

(3) 楼 的 直接 微 扰 方法 

将 微 扰 解 (2.39)-(2.40) 写成 如 下 形式 


A = 2AsechZe~*, (2.51) 
Z=2(X — 8), (2.52) 
e(z,t) = 2k(X -- £) + ô. (2.53) 


则 
X = ze^?d, (2.54) 
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1 = Ibe“, (2.55) 
k = ko, (2.56) 
+ koe ^e 
"ο. (257) 
2 L25-4« 2... J2,—4et 2 
σα κα es Ls (2.58) 
€ € 2 
当 E 一 0， 有 
X =r, 
€ = 如 τ Akot, 
δ = ôo — 4(k 1-2), (2.59) 
l= lọ. (2.60) 


此 时 ， 微 扰 解 退化 为 无 微 扰 NLSE 的 单 孤子 解 。 进 一 步 地 ， 比 较 不 同 微 扰 方 法 
得 到 的 E δ, X 等 参量 ， 不 难 发 现 本 论文 中 所 用 的 直接 微 扰 方法 得 到 的 近似 
解 中 的 各 个 参量 比 其 他 两 种 方法 得 到 的 具有 更 高 阶 的 修正 。 而 且 ， 利 用 该 直接 
微 扰 方法 ， 我 们 还 得 到 了 近似 解 的 一 阶 修正 。 所 以 ， 我 们 有 理由 相信 ， 利 用 楼 
的 直接 微 扰 方法 得 到 近似 解 有 更 高 的 准确 性 。 


(=) 直接 微 扰 方 法 的 可 靠 性 分 析 
上 一 节 中 ， 我 们 详细 叙述 了 楼 提出 的 直接 微 扰 方法 ， 并 从 与 其 他 微 扰 方 法 
得 到 的 近似 解 的 比较 中 发 现 ， 该 直接 法 得 到 的 解 具 有 更 好 的 准确 性 。 为 了 进 一 
步 从 量 上 来 确定 该 直接 微 扰 方 法 的 可 靠 性 。 这 一 节 中 ， 我 们 将 直接 微 扰 方法 运 
用 到 可 积 的 含 修正 项 的 NLSE 来 获得 它 的 近似 解 ， 借 助 Maple 工具 将 得 到 的 
近似 解 与 其 精确 解 进行 比较 和 分 析 ， 以 此 来 确定 直接 微 扰 方法 的 可 靠 性 701, 


考虑 如 下 的 含 修正 项 的 NLSEI7， 


it + iruz + uss + o|ul?u = 0. (2.61) 


C) wm 23- 


其 中 下 标 表示 对 时 间 变 量 t 和 空间 变量 z 的 求 导 。 在 非 线性 光学 意义 下 ， wu 
示 光 脉冲 的 缓 变 包 络 振幅 函数 ，z 和 + 分别 表示 脉冲 在 光纤 中 的 传输 距离 和 时 
间 ， 方程 中 y = 0k/dw = 1/v, (v ABBE), = B2k/Bw?， 由 于 光波 的 传 
播 常数 k 与 频率 w 的 依赖 关系 ， 不 同 频率 的 波 其 传播 速度 会 不 同 ， 由 此 将 产 
EGR, Hy, n 是 反映 光纤 的 色散 的 两 个 参量 。 方 程 中 第 四 项 描述 脉冲 的 非 
线性 效应 。 

1 精确 解 


采用 以 群 速度 v, 运动 的 参考 系 来 描述 方程 (2.61)， 即 
Z-z-tju 2z-ot T=t. (2.62) 
则 方程 (2.61) 可 简化 成 ; 
fur + puzz + o|ul?u = 0. (2.63) 


众所周知 ， 方 程 (2.63) 为 标准 NLSE。 它 有 如 下 形式 的 亮 孤 子 解 : 
u(Z,T) =% [TI tret Yr rasech [21(2 ~Zo+4kyT)}. (2.64) 


这 里 1, k, Zo, Oo 为 四 个 实 参数 ， 分 别 衡量 孤子 的 高 度 〈 宽 度 ) 、 速 度 、 初 始 
位 置 和 初始 相位 。 将 Z = z -mt，T =t ER, 185778 (2.61) 的 单 孤子 解 
为 (以 we 来 标记 精确 解 ): 


/2 
Ue = 2 = l ενα ome (E - Pol ech [o]( Ζ--γιί-- Zo + 4kyt)] 


2 i 
=% | ο fist m emp Mem Hlc ol (z — Zo + (4kw — )t). (3.65) 


此 时 ， 孤 子 的 高 度 (宽度 )， 初 始 位 置 和 初始 相位 仍然 由 参数 !，26，0o 确定 ， 
而 孤子 的 传播 速度 变 成 由 参数 k 和 γι 共同 表征 。 


2 近似 解 
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如 果 将 方程 (2.61) 中 色散 项 iyu 看 成 小 量 ， 我 们 就 可 以 用 直接 微 扩 方法 
对 其 进行 处 理 。 令 γι = ch， 其 中 < 是 表征 微 扰 强 弱 的 小 参量 ，/ 为 任意 常数 。 
利用 直接 微 扰 方 法 ， 使 得 方程 (2.61) 的 解 u 具有 如 下 形式 : 

u= HOU, re) = HD (uc, τ) + eu(C, 7) +--+). (266) 


其 中 a = a(z,t), 6 = B(z,t), Cm Cte, T= r(zte) 均 为 待定 的 实 函 
数 ， ς, T 有 如 下 性 质 ， 
(67) £3 (2,1). (2.67) 
XXE, Me 一 0 时 ，wo 退化 为 为 NLSE(2.63) 的 精确 解 。 
将 展开 式 (2.66) 代入 方程 (2.61)， 并 令 e 的 相同 次 智 的 系数 等 于 零 ， 得 到 
如 下 各 级 近似 方程 ; 
ior + "j2uoccC2 + e" ^o |uo[^us = 0. (2.68) 
ἱευχττι + peuk + 2e" oe|uo|"u, + e" oeudut + icaro 
--εβειιο + 1260, uo + ύγρεβ..ᾳ + χγρεασυαςς, + 2ίγρεβ-ιιος(, 
Γίεμυοςς, + Yotioc Ces + iuocki + "oor T4; + 2γοάρτςτιζ, = 0. (2.69) 


我 们 知道 ， 方 程 (2.68) 中 wo 是 不 显 含 e 的 ， 所 以 变量 C, 7 依赖 于 函数 α 
的 关系 式 为 : | 
σεν qup (2.70) 
并 且 uo 为 以 下 NLSE 
fuo, + Ύρυοςς + o|uo|^uo = 0 (2.71) 
的 精确 解 ， 是 方程 (2.71) 的 线性 化 方程 
iuir + Ύρύηςς + 20 Juo u Tow =0 (2.72) 
的 解 。 它 的 最 简单 的 两 个 形式 可 以 表示 成 ; 


Uy = Uor, Ur = Ue. (2.73) 


=) 垃 接 微 扰 方法 的 可 靠 性 分 析 ，25 - 


给 定 一 个 非 平凡 解 us, HT (2.69), (2.72) 和 (2.70) 对 于 任意 z, t 都 是 相 
容 的 ， 所 以 方程 (2.69) 中 其 他 项 的 系数 应 全 部 为 等， 即 


to, — βι + Oz, + ibzz = 0, ‘bz7z = 0, 


2εα;ς; + 2ieB,C. + deuce + Czz = 0. (2.74) 


解 之 得 
a=C, B= —pt/2, r=te®, C= ze. (2.75) 


其 中 C 为 任意 常数 。 因 此 ， 方 程 (2.61) 的 微 扰 解 可 表示 成 以 下 形式 (以 ua 标 
记 近 似 解 ): 
Ug = eO P e + eu). (2.76) 


3 ”可 靠 性 分 析 


图 2-1 给 出 了 不 同位 置 处 精确 解 的 强度 juu [P 和 近似 解 的 强度 juu? 随时 
间 的 演化 图 。 其 中 近似 解 中 wo 取 如 (2.64) 式 形式 ， 参 数 为 : (8) & = 0.5, l= 
06 C = 0.8, y = e€ = 0.01, Z% = % = 0, p = % = 1, o = 2 fl (b) 
k=0.5, L= 0.6, C = —04, yı =€=0.5, Ζο-θο-0, u234 21,0 —2. Æ 
传播 过 程 中 ， 两 列 孤 立波 分 别 保持 各 自 的 高 度 、 速 度 、 相 位 不 变 ， 这 正 是 孤立 
波 的 一 个 重要 特征 。 当 e 足够 小 时 ( 见 图 2-1(9)), FE (2.61) 的 近似 解 与 精确 
解 符合 得 相当 好 。 这 时 ， 直 接 微 扰 方 法 给 出 的 近似 解 有 效 地 描述 了 方程 (2.61) 
的 解 的 情况 。 不 过 引进 直接 微 扰 方 法 改变 了 孤 波 的 传播 速度 ， 此 时 精确 解 比 近 
似 解 传播 得 快 ， 随 着 传输 距离 z 的 增加 ， 近 似 解 开始 稍微 落后 于 精确 解 。z 越 
大 ， 近 似 解 偏 移 精 确 解 的 程度 也 越 明显 。 另 外 ， 引 进 直接 微 扰 方 法 增加 了 参数 
C， 近 似 解 的 强度 的 大 小 对 它 也 有 一 定 的 依赖 关系 ， 通 过 选取 合适 的 C 可 以 给 
出 与 精确 解 较 符合 的 近似 解 。 

^ e 值 增加 时 ( 见 图 2-1(b))， 两 列 孤 立波 的 传播 速度 都 减 小 (它们 的 具体 
形式 可 由 解 (2.65) 给 定 )， 其 中 近似 解 的 孤子 落后 于 精确 解 中 孤子 的 传播 ， 这 
就 使 得 在 z 逐渐 增加 的 过 程 中 近似 解 与 精确 解 越 来 越 不 匹配 。 此 时 ， 一 阶 修正 
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COO9 Exact solution 
—— Approximate solutior 


图 2-1: 不 同位 置 处 精确 解 和 近似 解 的 振幅 随时 间 的 演化 图 , 参数 分 别 为 
(àk = 05,1206 C = 0.8, = € = 0.01, % = Zo =0, p= =1, 0 =2. 
(b)k = 0.5, [20.6, C = —0.4, 7, =€=0.5, bo = 2 =0, 155 21,072. 


Fig.2-1: The time evolution of amplitudes of one-soliton solutions with different dis- 
tances. The parameters are (a)k = 0.5, 12 0.6, C = 0.8, yı = e = 0.01, % = Zo = 
0, =Y =l, o=2 (bl 206, C = —0.4, k = =€=0.5, bo = Zo = 0, p= y = 
1,022. 


(四 ) 微 扰 对 孤子 的 影响 - 27- 


u 相对 于 零 级 解 wo 已 经 不 再 是 小 量 ， 而 且 近 似 解 (2.76) 对 于 方程 (2.61) 也 不 
再 是 一 个 有 效 的 近似 解 。 


R 21: z = 10 处 相同 时 间 下 近似 解 的 最 大 强度 [us |? 与 精确 解 的 最 大 强度 uL]? 的 相 
对 误差 的 比较 图 (其 中 jue? = 1.44) 


X 2-1 给 出 了 z= 10 处 相同 时 间 下 取 不 同 的 ε 值 时 近似 解 的 最 大 
强度 luol? 与 精确 解 的 最 大 强度 |u|? 的 相对 误差 的 比较 。 其 中 参数 取 为 
k=05,1=06, Zo = 如 =0, k= = 1,0=2， 此 时 精确 解 的 最 大 强度 为 
lu? = 1.44。 表 中 [Al = ||υε|} — [us[?| 表示 精确 解 的 最 大 强度 与 近似 解 的 最 大 
强度 的 绝对 误差 。 从 表 中 我 们 观察 到 ， 当 e 小 于 0.1 时 ， 相 对 误差 均 小 于 5%。 
一 般 地 ， 我 们 认为 此 时 直接 微 扰 方法 可 以 很 好 地 描述 微 扰 方 程 的 情况 。 随 着 e 
的 增加 ， 两 者 之 间 的 误差 也 有 所 增加 。 

通过 以 上 的 简单 分 析 ， 我 们 发 现 当 参量 e 值 相 对 于 所 要 研究 的 问题 来 说 足 
够 小 时 直接 微 扰 方 法 可 以 很 好 地 描述 微 扰 方程 的 解 的 情况 ， 此 时 该 直接 微 扰 方 
法 对 于 我 们 研究 微 扰 NLSE 中 孤子 的 传播 有 很 好 的 可 靠 性 。 而 且 因为 微 扰 普遍 
存在 ， 所 以 我 们 所 采用 的 直接 微 扰 方法 对 研究 微 扰 对 NLSE 中 孤子 的 影响 有 广 
泛 的 意义 。 


(四 ) 微 扰 对 孤子 的 影响 
AX (2.39)-(2.40) 中 看 到 ， 由 于 受到 微 扰 项 —ieA 的 影响 ， 孤 子 的 振幅 随 
着 时 间 的 增加 指数 衰减 ， 而 且 相 位 也 发 生 了 改变 。 进 一 步 地 ， 令 (2.39) 式 中 


αἱ —4et 
sech| alk LM e 


2lxe^?* + Lao] = sechó, (2.77) 
并 将 φ EFE e Wt, Bn 


ϕ 55 2l[(—1 + 2et)(4kt + α) + αρ]. (2.78) 
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Ei, RIER, WSR v = 4k 并 没有 发 生 改 变 ， 不 过 波 数 p= 
21(-1--2ed) 不 再 是 常数 ， 而 是 随 着 时 间 t 的 变化 而 变化 。 


(a) 0) 


Nth 
η s 


ANNON 


2-2: (a)NLSE (118918 (2.2) 的 强度 演化 图 . (b)uo W (a) 中 相同 
形式 时 ， 微 扰 NLSE(2.25) 的 近似 解 (2.39) 的 模 方 |4| 演化 图 . 
(2.2)-(2.3) 式 中 的 参数 取 为 : 有 = 0.25, 1 =0.5, zo = do — 0. 

Fig.2-2: (a)Evolution of intensity profile of one-soliton solutions of NLSE. 
(b)Evolution of intensity profile |A|? of solution of perturbed NLSE with e = 0.06. 
The parameters in (2.2)-(2.3) are k = 0.25, 1 = 0.5, zo = ġo = 0. 


2-2(a) 中 ， 我 们 给 出 了 NLSE 的 单 孤 立波 (2.2) 的 强度 的 演化 ， 图 2- 
2(b) 则 为 取 up 与 2-2(a) 图 中 相同 形式 时 ， 微 扰 NLSE 的 近似 解 (2.39) 的 模 方 
|Al? 的 演化 图 。 方 程 (2.2) 和 (2.3) PAAR k= 0.25, 1 = 0.5, zo = fo = 0 W 
及 (2.25) RP e= 0.06。 从 2-2(a) 图 我 们 知道 ， 当 没有 微 扰 作用 时 ， 孤 立波 沿 
着 > 轴 稳 定 地 传输 ， 它 的 形状 、 传 播 方 向 和 速度 都 没有 发 生变 化 。 但 是 由 于 存 
在 损耗 ( 即 微 扰 项 -ic4)， 珀 子 的 形状 发 生 了 明显 的 变化 : Ht, MERES 
持续 减少 使 得 孤子 的 高 度 随 着 时 间 的 增加 作 指 数 衰减 。 当 + 一 oo， 孤子 振幅 也 
逐渐 减 小 并 趋向 0。 而 且 很 明显 地 ， 孤 子 的 传播 方向 也 发 生 了 改变 。 

2-3(a) 表示 NLSE 的 双 孤 立波 的 强度 演化 图 。 取 wo 与 2-3(a) 图 中 相 
同形 式 时 ， 微 扰 NLSE 的 近似 解 的 模 方 Al? 随时 间 的 演化 如 图 2-3(b), X 
f£ (2.25) PER e = 0.01。 当 没有 微 扰 作用 时 ， 两 孤立 波 沿 着 x 轴 稳 定 地 传输 ， 
它 的 形状 、 传 播 方向 和 速度 都 不 发 生变 化 。 但 是 从 图 2-3(b) 我 们 观察 到 ， 微 扰 
—ieA 使 孤子 的 形状 发 生 了 明显 的 变化 ， 此 时 ， 孤 波 的 振幅 也 是 随 着 时 间 的 增 
加 而 指数 衰减 并 逐渐 趋 于 0。 但 是 ， 微 扰 并 不 会 改变 两 孤子 间 的 能 重 分 布 ， 即 


(a) 0) 


Bl 2-3: (a)NLSE 的 双 孤 波 的 强度 演化 图 ，(b) 取 wo 与 (e) 图 中 相同 形式 且 
<= 0.01 M, (Kit NLSE(2.25) 的 近似 解 的 模 方 LAU? 的 演化 图 . 
(2.11)-(2.16) 式 中 取 kı = 0.5 +0.2i, kp = 0.5 — 0.5i, 1 = nf — 0. 
Fig.2-3: (a)Evolution of intensity profile of two-soliton solutions of NLSE. 
(b)Evolution of intensity profile |A|? of solution of perturbed NLSE with ¢ = 0.01. 
The parameter in (2.11)-(2.16) are kı = 0.5 + 0.2i, kz = 0.5 — 0.57, no) = η) =Q. 
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地 分 布 在 两 分 量 上 而 不 转移 或 交换 。 所 以 说 ， 碰 撞 的 性 质 不 会 因为 微 扰 的 影响 
而 发 生变 化 。 


(五 ) ^5 


利用 直接 微 扰 方法 处 理 NLSE 的 微 扰 问 题 的 过 程 中 得 到 的 近似 解 中 的 零 阶 
解 可 以 取 微 扰 NLSE 的 精确 解 ， 而 且 是 任意 一 个 ， 而 一 阶 修正 则 可 以 为 无 微 
扰 方 程 的 任意 一 个 对 称 。 众 所 周知 ，NLSE 的 精确 解 和 对 称 均 有 无 穷 多 个 。 因 
此 ， 采 用 楼 的 直接 微 扰 方法 得 到 的 近似 解 明显 比 利 用 其 他 的 微 扰 方法 得 到 的 解 
要 丰富 得 多 。 这 里 ， 相 同 的 等 阶 解 与 不 同 的 一 阶 修正 可 以 组 成 不 同 的 近似 解 ; 
或 者 ,不同 的 零 阶 解 和 相同 的 对 称 也 可 以 组 成 不 同 的 近似 解 。 

2-4 给 出 了 在 不 同时 间 下 ， 取 相同 的 等 阶 解 与 两 个 不 同 的 一 阶 修正 
组 成 的 两 个 近似 解 的 振幅 之 闻 的 比较 ， 其 中 (a) 为 单 孤子 ; (ο) 为 双 孤 子 。 
图 中 用 圆圈 表示 一 阶 修正 取 u = woe 的 情形 ， 实 线 表示 u = uon. w 则 分 
别 由 (2.2) 式 (MMF) 和 (2.11)-(2.16) X. (MMF) 来 确定 。 式 中 参数 分 别 
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Rk =0.25,1=0.5, zo —0, €= 0.2 Mk, =054+0.2i, ka = 0.5 —0.5i, 
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图 2-4: 相同 的 等 阶 解 不 同一 阶 修正 情况 下 两 个 不 同 近似 解 的 模 方 的 比较 图 . 
图 中 贺 轿 表示 一 阶 修正 取 wi = woe， 实 线 表示 u = ug. (a) MMF; (b) MHF. 
Fig.2-4: The time evolution of amplitudes of solutions with the same zero-order 


solution and different first-order solutions. (a) one-soliton; (b)two-soliton. 


"D =n = 0，e = 0.2。 我 们 知道 ， 本 身 一 阶 修正 并 不 是 对 解 起 主要 作用 的 
部 分 ， 它 的 改变 只 是 很 微小 的 。 反 映 在 图 2-4 中 的 结果 就 是 ， 取 不 同 的 一 阶 修 
正 时 ， 两 个 近似 解 的 振幅 几乎 完全 重合 ， 虽 然 此 时 的 e = 0.2( 已 经 不 算是 小 量 
了 )。 
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上 一 章 我 们 以 (1+1) 维 的 微 扰 NLSE 为 例 ， 详 细 叙 述 了 楼 提出 的 直接 微 扰 
方法 ， 并 继而 将 采用 直接 微 扰 方法 获得 的 近似 解 与 用 其 他 微 扰 方 法 得 到 的 解 进 
行 了 简单 的 比较 。 从 中 不 难 发 现 ， 楼 的 直接 微 扰 方 法 有 较 大 的 优越 性 ，(1) 不 
依赖 逆 散 射 方 法 ， 也 就 是 在 这 种 直接 微 扰 方 法 中 ， 不 要 求 无 微 扰 的 NLSE 是 可 
以 用 道 散 射 方法 精确 求解 的 ，(2) 思路 直接 ， 容 易 理 解 和 接受 ; (3) 利用 楼 的 直 
接 微 拢 方法 能 同时 得 到 各 种 不 同类 型 的 微 扰 解 ， 包 括 单 孤 子 型 ， 多 孤子 型 ， 周 
期 型 解 等 等 。 在 解决 了 (1+1) 维 NLSE 的 微 扰 问 题 后 ， 我 们 自然 想 问 ， 该 直接 
微 扰 方法 能 不 能 用 来 解决 其 他 的 NLSE 的 微 扰 问题 ， 如 高 维 的 NLSE, BAE 
Rete CIs ARS Ss? 

据 我 们 所 知 ， 到 目前 为 止 大 多 数 的 微 扰 方 法 均 倾 向 于 处 理 单 场 量 的 非 线性 
演化 方程 的 微 扰 问题 ， 很 少见 到 有 处 理 耦 合 非 线性 演化 方程 组 的 微 扰 问 题 的 ， 
而 鉴于 它们 在 各 种 不 同 的 物理 领域 的 广泛 应 用 ， 它 们 已 经 成 为 人 们 广泛 关注 的 
课题 。 另 外 在 实际 问题 中 还 存在 着 大 量 的 二 维 甚至 更 高 维 的 非 线性 演化 方程 ， 
关于 它们 的 微 扰 问 题 ， 迄 今 也 仍然 难于 求 得 其 解析 形式 的 解 ， 而 只 能 通过 计算 
机 去 求 数值 解 或 进行 模拟 。 

所 以 在 这 一 章 中 ， 我 们 将 证 实用 直接 微 扰 方法 也 能 解决 耦合 非 线性 薛 定 订 
方程 组 以 及 (2+1) 维 NLSE 的 微 扰 问题 178, 791, 


(一 ) BARER HEA 
HARRESA (CNLSE) 
iut + Οὐχ + (antul? + aazlu|2)v 一 0, 


ive + ὀχυεα 十 (aazlo + a; |u[?)v = 0. (3.1) 


作为 描述 非 线 性 现象 的 基本 模型 之 一 ， 它 出 现在 许多 不 同 的 物理 领域 中 ， 如 非 
线性 光学 80, 811， 生 物 物理 8 为 ， 凝 聚 态 物理 83, 84) 等 ， 因 此 也 越 来 越 受到 人 
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们 的 关注 。 尤 其 地 ， 在 非 线 性 光学 中 ，CNLSE 用 于 描述 双 折 射 光纤 中 不 同 模 
式 的 光纤 的 相互 作用 ， 其 中 u, v 表示 两 相互 作用 的 光纤 模式 的 缓 变 振幅 ，c1 
和 c 为 色散 系数 。 这 个 方程 既 包 括 了 自 相 位 调制 (SPM) juu 和 |u?v, ΧΕΙ, 
含 了 交叉 相位 调制 (XPM) lvu 和 |wl2v。 同 相位 调制 是 形成 孤子 的 先决 条 件 ， 
而 交叉 相位 调制 则 引起 不 同 模式 (包括 极 化 ) 或 不 同 信道 之 间 的 耦合 。 当 两 束 或 
更 多 束 光波 同时 在 光纤 中 传输 时 ， 它 们 将 通过 光纤 中 的 非 线性 发 生 相 互 作 用 。 
通常 ， 这 样 的 相互 作用 在 适当 的 条 件 下 通过 不 同 的 非 线性 现象 会 产生 新 的 波 ， 
如 受 激 喇 曼 散 射 和 受 激 布 里 渊 散 射 、 谐 波 产生 以 及 四 波 混 频 ， 但 是 ， 相 同 的 非 
线性 通过 上 述 交 叉 相 位 调制 现象 也 可 在 入 射 波 之 间 产 生 耦 合 。 

RE 20 世纪 70 年 代 ，Manakovl82 首次 利用 IST 研究 了 CNLSE (3.1) 
的 可 积 性 ， 他 发 现 当 a = o = 1 an = an = Og = On 时 (满足 该 条 件 
的 CNLSE(3.1) 后 来 也 因此 被 称 为 Manakov 方程 )， 方 程 (3.1) 是 可 积 的 。 后 
X, Zakharov 和 Schulman?! 关于 方程 (3.1) 的 可 积 性 的 研究 对 Manakov 的 
结果 进行 了 补充 ， 他 们 认为 当 且 仅 当 


Cr = Ca, 041 三 022 = An = an (3.2) 
或 者 
€; — —09, Q = 02? = 一 Q12 = 一 Q21 (3.3) 
时 ， 方 程 组 (3.1) 才 是 可 积 的 ， 即 它们 可 以 用 IST 进行 求解 ， 除 此 之 外 的 
其 他 的 情况 都 是 不 可 积 的 。 同 样 的 可 积 条 件 也 曾 被 Radhakrishnan 等 人 利 
用 Painlevé 分 析 法 报道 过 (841, 
文献 [82] 中 ，Manakov 还 进一步 给 出 了 CNLSE(3.1) 的 单 孤立 子 解 : 


u = cos(P)e peer 6-96 gocher — 2pgt + 70), 
v = sin(0)e /2petes t") od sech(pz — 2pgt + γρ). (3.4) 


其 中 p, q, 0, δ, Yo, Go 都 是 任意 实 常数 。 从 解 (3.4) F, RUT MM I 
子 的 振幅 V2p， 速 度 2g， 初 始 的 中 心 位 置 -7o/p 等 参数 。 

继 Manakov 的 开拓 性 的 工作 后 ， 越 来 越 多 的 物理 学 家 为 CNLSE 的 解 的 
研究 投入 了 大 量 的 精力 ， 特 别 是 它 的 孤立 子 解 。 如 Radhakrishnan 利用 Hirota 
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方法 给 出 了 Manakov 方程 的 单 孤子 解 1, 


Qe 
(lrentueR] 


pen 


ig + επι τη Β΄ 


u 


其 中 
,a Mo + (6?) 
(ki + ky)?’ 


η) = y(t + ἐκ) +. 
双 孤子 解 (86. 871, 


04€"! + age? + emtnitmtS: 十 οὖν η Γη 53 


u= non, aa LESS a ee 
1+ emtnith: + emtnptSo + enitmtss 十 emn +R 4 emtnitmtny+Rs’ 
Bien + Boe™ 十 emtnitmts, 4 emtm+nj+S; 
m LM——————— M ————M——————— συ — —— 
1+ οτι ΓΙ ΓΙ + οι ΓΠΣ 50 十 enim S6 十 erani Ra 十 em-eni-emeni Rat 


Qs = cos(0, e», bs = βἰη(θ,)ετύδ», 


Ne = aE + ik2T + nao, ks = p, + igs, 


Nso = Yeo + 16,0; eg = rar 
2 
Rm _ M41 εἴθ 一 K22 
hr ka + E 
k-k 
$i 1 2 = 
e [EST ES F ΚΠ (ky n kx) (aika 03K11), 
" k-k 
e? = Eres πο NNUS 一 - 
(ky 十 kt) (kt ef k2) (81521 fata), 
δα _ ko — ky ην 
€ (he + (ki +) + ks) (ky ip sts 022), 
ko — k 
S 2 M . - 
[ESTIS xiu Bike), 
Rs lk — kj? 


x: (ht (ky +S), + Bye — Ki2K21), 


. 33 . 


(3.5) 


(3.6) 


(3.7) 
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Kij — 2(αι + o7) . 
这 里 Ps, ds; bs, δε, σου Ceo (s = 1,2) 均 为 任意 实 常数 。 


(3.8) 


1 近似 解 
为 简单 和 方便 起 见 ， 在 本 节 中 我 们 只 研究 当 er = o = an = an = an = 
025 = 1 时 的 CNLSE(3.1) 的 微 扰 问题 ， 即 
fu, + Use + (Jul + |vP)u = —ieu, 
fu, + Όχα + (lvl? + luf?) = —iev. (3.9) 
在 非 线 性 光学 中 ， 微 扰 项 ieu 和 dev 表示 线性 损耗 ， 其 中 ε 为 衡量 微 扰 强度 的 
小 参量 。 当 e = 0 时 ， 方 程 组 (3.9) 将 退化 为 标准 的 CNLSE(3.1). 
利用 楼 的 直接 微 扰 方法 ， 将 方程 组 (3.9) 中 的 和 vo 展开 成 ; 
u = etu (E, r, ε) = eft) lun (€, τὴ + eu(é, 7) + Ο(ε2)], 
v= eeu (E, r, e) = ee volt, τ) + ev; (£, 7) + O(e)]. (3.10) 
这 里 ， 
a=a(z,t), b = b(z,t), cz c(z,t), 
dz d(z,t), € = €(z,t,e), n=n(z,t,€) (3.11) 
均 为 待定 的 实 函 数 ， 并 且 变 量 {é&,7} 满足 以 下 关系 式 : 
{Εντ} 5 {r,t}. (3.12) 
因此 ， 当 ε-θ 时 ， up 和 Uo 将 退化 为 无 微 扰 的 CNLSE 的 精确 解 。 
将 展开 式 (9.10) 代入 微 扰 的 CNLSE (3.9)， 得 到 
itor Te + Uinger” + e° uo uo + e° lvo uo + elitt: + ueta” 
+e" (uuo + 2|uo|?u;) + e (u |vo|? + υφυιυξ + uovovt)] 


tielas + ἴδε)ωο + e(azz + ibzz)uo + ieuo + i£ruoe + ste 
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θείας + ibz)Ecting + 2Εετε ρε; + Tao = Ο(ε”), (3.18) 
ivorte + VogeEx” + e lvo vo + ezeejuo|2oo + ε[έυιττι + Viggen” 
"e" (vivo? + 2/υο]”νι) + e% (vuo? + vous us + vououi;)] 
+ie(c: + íd;)ug + € (Con + idu, )vg + devo + evoe + ξεσυηε 
+2€(Cy + id; )EaUog + 2eTaVoer + Τοσυοτ = O(€). (3.14) 


由 于 ér 有 如 (3.12) 的 性 质 ， 所 以 &, το 和 6. 都 是 e 的 一 阶 项 。 令 方程 
组 (3.13) 和 (3.14) Pe 的 相同 次 容 的 系数 均等 于 零 ， 得 到 如 下 各 级 近似 的 方程 
组 ; 

fug, Tt + Uoeeéz? + €" ^ uol uo 十 ezeeluol2uo = 0, 


Wort + υρεεξε” + e lvo] vo + €* uol? vo =0. (3.15) 


ε[ίυηττι + ugein + ο ^ (utug? + 2|uo[^u;) + e (ulvo? + uous 
-Fugugt; )) + ἱε(αι + ibi)uo + (a. + ib, )ug + euo + irung 
Ἑξεείίος + 2ε(α» + δε) δεις + 26, Totlogr + Trztor = 0, (3.16) 
ε[ίνχντι + Veeés’ + e" (vivo? + 2|vo]?vi) + e" (vi |uo[? + vou, ug 
-Fuuou; )] + ἑε(οι + id;)vo + ε(οες + ἱάεα)υο + ievo + ἑξγυας 
TÉzz Uoc + 2Ε(Ος + id; )ExVoe + 26, T, vor + Tz Uo, = 0. (3.17) 


我 们 知道 ， 方 程 组 (3.15) 中 的 场 量 w 和 w RES e， 所 以 变量 ér 应 满 
RRA: 
ἔν =e* me Al n=e - gr, (3.18) 
解 之 得 c = co KH, ABA (3.15) 化 为 无 微 扰 的 CNLSE( 不 过 此 时 ， 方 程 组 
以 6,7 为 变量 ): 
fuo. + woee + (luol? + lvol?)uo = 0, 


Wor + Voge + (Ivo? t |uo|?)ug = 0. (3.19) 
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从 展开 式 (3.10) 我 们 知道 ， M v 正 是 CNLSE (3.19) 的 对 称 ， 即 它们 是 方 
程 组 (3.19) 的 线性 化 方程 组 
iur + Urge + uiuo? + luo u + ulvo]? + ugvi vo + uovov? = 0, 
ivir + Vige + UJ vo? + 2|vo[^vi + vi |uo[? + vou; us + vougut = 0 (3.20) 
的 解 。 它 的 最 简单 的 两 个 解 可 以 表示 成 831， 
Ui = Ue, U1 = Voge 和 uy = Uor, Uy = Vor. (3.21) 


给 定 一 个 非 平凡 解 up 和 veo， 方程 组 (3.16), (3.17), (3.18) 和 (3.20) 对 
于 任意 ct 都 是 相 容 的 ， 所 以 方程 组 (3.16) 和 (3.17) 中 其 他 项 的 系数 应 全 部 
为 等 ， 即 


ber --αι +1=0, au ~b 0, (3.22) 
Erz + 2eaz€,=0, & + 2b, Ξ 0, 7, — 0, (3.23) 
dj +G¢+1=0, c,—d, =0, (3.24) 
Ezz + 2εΟρξα — 0, ἕι ἠ- 2εάςξ, — 0, τες, (3.25) 


解 方程 组 (3.22)-(3.25) 和 (3.18), ABBE a, b, ο, d, & τ 的 表达 式 如 下 : 


a=c=-—2t, bes +o 


s H 2 5 —2%t NN 1- e^t 
d= 52 +o, =r“, r= FE (3.26) 
其 中 οι, ο 均 为 任意 实 常 数 。 
这 样 ， 微 扰 的 CNLSE(3.9) 的 近似 解 可 以 表示 成 : 
u- e-i) (uo + eu), 
v = ett) (to + ev). (3.27) 


其 中 我 们 取 οἱ = ο) = 0ο 


2 微 扰 对 孤子 的 影响 
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到 目前 为 止 人 们 已 经 找到 了 无 微 扰 的 CNLSE(3.19) 的 各 种 各 样 的 解 的 形 
式 ， 尤 其 是 孤立 子 解 ， 包 括 单 孤立 子 解 和 多 孤立 子 解 等 。 而 近似 解 (3.27) 中 的 
一 级 修正 u Mu 则 可 以 由 (3.21) 式 确 定 。 在 此 基础 上 ， 本 节 中 我 们 将 结合 微 
H CNLSE 在 非 线性 光学 中 的 应 用 来 讨论 损耗 项 ieu 和 —iev 对 孤子 的 各 个 参 
数 产生 的 影响 。 

从 第 一 章 中 对 (1151) 维 NLSE 的 近似 解 的 分 析 中 我 们 已 经 知道 ， 解 (3.27) 
PSOE uo 和 vo 也 可 以 是 无 微 扰 的 CNLSE(3.19) 的 任意 一 个 精确 解 ， 而 
一 阶 修正 u 和 vw 则 可 以 是 (3.19) 式 的 任意 一 个 对 称 。 图 3-1 描述 了 CNLSE 
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图 31: 参数 取 p=0.5, «ΞΟΞΣΞΙ,ΘΞ Ξ, (020 时， 
单 孤立 波 (3.5) 的 强度 (a)|uo|? 和 (b)|vo[? 随时 间 的 演化 图 ， 
Fig.3-1: Evolution of intensity profile (a)|uo|? and (b)|vo|? of one-soliton solutions of 


CNLSE. The parameters are fixed at p = 0.5,q= % = b= aes 5, 40 = 0. 
9 


图 5-2: SBME uo 和 vo 取 与 上 图 相同 形式 时 ， 
微 扰 CNLSE 的 解 (3.27) 的 模 方 (cjjw|> 和 (d)|v|? 随时 间 的 演化 图 . 
Fig.3-2: Evolution of intensity profile (c)|u|? and (d)|v|? of solution of perturbed 
CNLSE with e = 0.08. The up and vp in (3.27) are the same as those in Fig.3-1. 
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的 单 孤立 波 的 强度 (auol? 和 (b)luo|? 随时 间 的 演化 情况 。 此 时 ， 两 列 孤立 波 
分 别 沿 着 各 自 的 α 轴 作 稳定 传播 而 不 改变 波 的 形状 、 传 播 速 度 和 方向 ， 它 们 的 
高 度 、 传 播 速度 、 初 始 位 置 和 初始 相位 分 别 由 p、 σα. mw Co 等 参数 决定 。 当 
有 微 扰 作用 时 ( 见 图 3-2)， 孤 立波 的 强度 随时 间 的 演化 与 (1-1) 维 NLSE 的 情 
况 非常 类 似 。 此 时 ， 孤 子 的 振幅 和 能 量 都 随时 间 的 增加 而 指数 衰减 直至 零 。 微 
16535 e 值 越 大 ， 衰 减 得 也 就 越 快 。 另 外 ， 由 于 微 扰 -iew 和 一 iev HH, M 
立 子 的 相位 和 波 数 等 也 都 发 生 了 相应 的 改变 。 它们 的 具体 形式 可 以 通过 近似 


(a) ο 


if P 
ly 


图 3-3: 非 弹 性 碰 擅 的 双 孤 立波 的 强度 (a)|uo|? 和 (b)jvo[? 的 演化 图 . 
参数 取 ky = 0.5-F0.5i, ko = 0.55 — 0.25i, ay = az = By = 1, f = 0, mo = nzo = 0. 
Fig.3-3: Evolution of intensity profile (a)|uo|? and (b)|vo|? of two-soliton solutions of 
CNLSE. The parameters are fixed at kj = 0.5 + 0.57, ko = 0.55 — 0.253, 


Qj = 02 = (i = 1,  — 0, mo = Nao = 0. 
(o (4) 


Β 3-4: 零 阶 和 解 uo 和 vo 取 与 上 图 相同 形式 时 ， 
微 扰 CNLSE 的 解 (3.27) 的 模 方 (c)juj? 和 (d)|v[? 的 演化 图 . 
Fig.3-4: Evolution of intensity profile (c)|u|? and (d)|v|? of solution of perturbed 
CNLSE with e = 0.02, and uo and vo in (3.27) are the same as those in Fig.3-3. 
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解 (3.27) 详细 地 表示 出 来 ， 有 兴趣 的 读者 可 参见 第 二 章 。 但 是 信 得 指出 的 是 ， 
孤立 波 始终 以 速度 Vip 传播 而 不 发 生变 化 。 

3-3 给 出 非 弹性 碰 擅 的 双 孤 波 的 强度 (a) |uo[? 和 (b)vo[? 随时 间 的 演化 
图 。 在 初始 情况 下 ， 两 孤立 子 沿 着 坐标 轴 稳 定 地 传播 ， 它 们 的 能 量 也 分 别 分 布 
在 两 个 孤立 子 上 而 没有 发 生 任 何 变化 。 但 是 由 于 交叉 调制 项 |uj2v 和 |v|?w 的 作 
用 ， 两 孤子 逐渐 靠近 并 发 生 碰撞 ， 此 时 两 孤子 上 的 能 量 不 再 分 裂 而 全 都 集中 在 
了 一 个 孤立 子 上 。 经 过 简单 的 计算 后 将 发 现 ， 孤 子 的 能 量 在 整个 碰撞 过 程 中 


图 3-5: SAW ky =1+i, ke Ξ-1 -- 0.5%, αι = 1 = az = fo = 1.5, mo = no = 0 Ff, 
弹性 磁 擅 的 双 孤 立波 的 强度 (a) |uo[? 和 (b)|vo[? 的 演化 图 . 
Fig.3-5: Evolution of intensity profile (a)|uo|? and (b)|vo|? of two-soliton solutions of 
CNLSE. The parameters are fixed at kj = 1 +i, kp = 1 — 0.53, 
αι --βι = 02 = =1, mo η. 
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图 5-6: 取 与 上 图 相同 形式 的 wo 和 vo 时 ， 
微 扰 CNLSE 的 解 的 模 方 (c)lwj? 和 (ἀ)|υ|2 随时 间 的 演化 图 . 
Fig.3-6: Evolution of intensity profile (c)|u|? and (d)lv|? of solution of perturbed 
ONLSE with e = 0.08. The ug and vo in (3.27) are the same as those Fig.3-5. 
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始终 守恒 。 这 一 性 质 在 文献 [86] 中 有 详细 描述 。 不 过 无 论 在 碰撞 前 还 是 碰撞 
后 ， 孤 立波 都 沿 着 α 轴 传 播 而 保持 各 自 的 高 度 和 速度 不 变 。 当 存在 线性 损耗 项 
时 ( 见 图 3-4)， 两 孤立 波 的 强度 随 着 时 间 的 演化 而 不 断 衰减 直至 零 ，e 值 越 大 ， 
衰减 得 也 越 快 。 


图 3-5 和 图 3-6 分 别 表示 CNLSE 的 弹性 碰撞 的 双 孤 立波 的 强度 演化 图 和 
微 扰 CNLSE 的 两 孤子 解 的 模 方 随时 间 的 演化 图 。 当 没有 征 扰 作用 时 ， 无 论 是 
碰撞 前 还 是 碰撞 后 两 孤立 波 随 着 时 间 的 增加 而 稳定 传播 。 与 非 弹 性 碰撞 情况 不 
同 的 是 ， 此 时 ， 珀 子 能 重 在 磁 撞 后 不 是 集中 在 一 个 孤子 上 ， 而 是 重新 又 分 布 在 
两 个 孤子 上 ， 且 各 自 的 能 量 保持 与 碰撞 前 的 一 致 而 不 发 生 任何 转移 和 交换 。 当 
有 微 扰 项 -ieu 和 —iev 的 作用 时 ， 两 孤子 的 强度 都 将 随 着 时 间 的 增加 而 持续 误 
减 直至 等 。 而 且 此 时 两 孤立 子 的 能 量 只 是 各 自 衰 减 ， 原 来 均匀 分 布 在 两 孤子 上 
的 能 量 不 会 因为 有 微 扰 影 响 而 重新 分 布 ， 也 就 是 说 ， 孤 立 子 的 碰撞 性 质 不 会 因 
为 微 扰 的 作用 而 发 生变 化 。 


为 了 进一步 检验 直接 微 扰 方法 的 有 效 性 ， 我 们 采用 分 步 傅立叶 方法 对 微 
δ CNLSE(3.1) 进行 了 数值 模拟 。 结 果 如 图 37 所 示 。 图 中 实 线 表示 利用 直接 
微 扰 方法 得 到 的 近似 的 解析 解 ， 虚线 为 数值 结果 。 图 3-7 中 (a) 表示 取 不 同 ε 
值 时 ， 微 扰 CNLSE 的 单 孤 子 解 的 振幅 随时 间 的 演化 图 。 近 似 解 (3.27) 中 的 堆 
阶 解 wo 由 (3.4) 式 给 定 ， 其 中 参数 为 p= 4-05, θ-π/4, ὃ--ηο-- ῷ--θ, 
且 取 一 阶 修正 v; = we。 从 图 中 可 以 观察 到 ， 当 e = 0.01， 解 析 解 与 数值 解 几 
FREES. ME e 值 的 增加 ， 两 者 之 间 的 差距 也 逐渐 增加 。 不 过 在 相当 大 的 
e 值 范围 内 ， 解 析 解 与 数值 解 易 合 得 相当 好 。 例 如 当 e = 0.1， 解 析 解 与 数值 解 
之 间 的 绝对 误差 仅 为 0.0012。 


图 3-7(b) 描述 双 孤 子 解 的 振幅 随时 间 演 化 的 情况 。 为 简单 起 见 ， 这 里 我 
们 选择 以 等 幅 的 两 孤子 为 例 来 说 明 情 况 ， 此 时 ， 能 量 平 均 地 分 布 在 两 孤子 
上 ， 所 以 它们 的 振幅 可 以 由 共同 的 一 条 曲线 完全 确定 。 零 级 解 w 由 (3.7)- 
(3.8) 式 给 定 ， 并 取 参 数 为 =1+i k = 1-054, αι = βι = mo = 1, 
o = f = V2/2, no = 2 以 及 一 阶 修正 u = uge WENA, RUT AE 
的 振幅 十 分 吻合 ，e 值 越 小 ， 两 者 符合 得 越 好 。 然 而 从 图 中 也 不 难 发 现 ， 解 析 
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解 和 数值 解 的 碰撞 时 刻 并 不 相同 ， 后 者 的 碰撞 时 刻 要 超前 于 前 者 的 。 特 别 地 ， 
在 碰撞 时 刻 解析 解 和 数值 解 的 振幅 也 有 明显 的 差别 : 随 着 e 值 的 增加 ， 解 析 解 
在 碰撞 时 刻 的 振幅 并 没有 发 生 明 显 的 变化 。 而 数值 结果 却 显示 ， 解 的 振幅 随 
18 ε 的 增加 而 明显 减 小 。 
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ΕΙ 3.7: 取 不 同 的 e 值 时 ， 孤 子 的 振幅 随时 间 的 演化 图 .(&) 单 孤子 (b) AMF. 
实 线 和 虚线 分 别 表示 解析 解 和 数值 解 的 振幅 . 
Fig.3-7: The time evolution of amplitudes of solutions with different € values. 
(a)one-soliton; (b)two-soliton. The solid lines are from analytical calculations, and the 


dashed ones are from numerical simulations. 
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Xd EAS GT) A MERER EIR PER 
现 ， 楼 的 直接 微 扰 方法 也 可 用 于 处 理 多 分 量 耦 合 的 非 线性 薛 定 廖 方程 的 微 扰 问 
题 。 

考虑 微 扰 的 多 分 量 耦 合 非 线性 薛 定 户 方 程 : 


iEms + CmEmtt + 20mm|Em|"Em +2 Y Amil Ei Em = —ieEm. 
(im) 

Ami = Qim, (m,l 2 1,2,--- N) (3.28) 
该 方程 在 实际 物理 问题 中 已 经 得 到 广泛 应 用 【79。 如 在 非 线 性 光学 中 ， 方 
程 (3.28) P Επι 为 光纤 的 第 πι 个 模式 的 慢 变 化 振幅 ，z, t 分 别 为 无 量 纲 化 后 
波 的 传播 距离 和 时 间 ， 参 数 ο 表示 色散 系数 ，amwm 为 波 包 的 自 相 位 调制 系 
数 ，a 描述 交叉 相位 调制 系数 。 微 扰 项 —ieE, 表示 光纤 中 的 线性 损耗 ， 其 中 
e 为 描述 微 扰 强度 的 小 参量 。 当 e = 0， 方 程 (3.28) 退化 为 无 微 扰 的 多 分 量 耦 合 
方程 组 


iE + CmEmtt + 20mm|Em|"Em + 2 )3 Omt| Ei] Επι = 0. (3.29) 
(πι 
一 般 地 ， 多 分 量 CNLSE(3.29) 是 不 可 积 的 ， 即 它 不 能 由 IST 精确 求解 。 
不 过 在 文献 [84] 中 作者 采用 Painlevé 分 析 法 对 方程 (3.29) 进行 可 积 性 研究 。 
他 们 发 现 ， 当 (3.29) 中 的 各 参量 满足 以 下 条 件 时 ， 该 方程 是 可 积 的 ; 
(1). c=o (§ =1,2...N) Qj, =Q. (,52 1,2... N) 
(2. -αξοιξο, (ὑπέἠ, j21,2...N) 
04, = (-1)55(—1)5*o. (1,8 — 1,2... Ν) 
K^ i-212...N (& i 只 能 取 一 次 ) 
(3. Cy = Cy = c; = 6, (hia j, j— 1,2... N) 
αι = [(—1)* (1) (-1)»»(-1)*«]a. (Ls 2 1,2... N) 
其 中 i, i2=1,2...N (i «i) (每 组 (i, ij) 只 能 取 一 次 ) 


(=) SP ERATAER EER - 43 - 


(n. 一 ca = Gig =...= Cy = CF = C, 
G Ai Ato x... Finn, J=1,2...N) 
04, = [(—1) (-1)5 ... (21): (-1)5-1]o. (ls = 1,2... N) 
其 中 i, izn .. 1 Ξ1,2...Ν (q«i «...«t) 
(每 组 (i, do, ...,in_1) 只 能 取 一 次 ) 
(nt1) | —6, 2 —64 =... = —6, = Cj = C, 
(σσ... σε, f=1,2...N) 
04, = [(--1) δε (--1) δν ... (—1) im (1) Ja, (ls = 1,2... N) 
其 中 i, to, ...νίη-1,9...Ν (ü«i«...«i) 
(每 组 (i, do, ..., i.) 只 能 取 一 次 ) 
其 中 若 N ABM, n=N/2; 若 N 为 奇数 ， 则 m = (Ν —1)/2 dy 满足 
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1 — 
利用 直接 微 扰 方法 ， 将 Em RAR: 
Ep = emtia) (pn (E, τὴ egn(€ A+++). (moL2.-N)  (&X) 
其 中 
am = Am(Z,t), bm = bs (zt), € = E(z,t,€), 7 =T(z,t,€) 


均 为 待定 的 实 函 数 。 
利用 直接 微 扰 方 法 ， 得 到 相 容 性 条 件 : 
és ze -- ε΄, τι e? = e, 


ία, — bme + Cman + icmbmit + i = 0, &n =0, 


Wm EAT, + 2iCmeDmeT + ὅτε + CmTu = 0. (3.31) 
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经 过 简单 的 运算 ， 不 难 求 得 : 
a=—22, bm =t?/2em, 


€=(1—e7*)/4e, T= te, (3.32) 
ΤΕ, JEA (3.28) 的 近似 解 可 表示 成 : 


Em = ο ο” /12om(pn + egm). (3.33) 
其 中 pm 为 无 微 扰 CNLSE 
N 
iE mg + cmEmrr + 20mm|Em|?Em +2 ^ am| ΕΙ} Em = 0 (3.34) 
{=1 
(žm) 


的 任意 一 个 精确 解 ，gm 是 它 的 任意 一 个 对 称 。 


2 微 扰 对 孤子 的 影响 


根据 文献 [86] 知 ， 当 cm =c=1, aim =a (l,m = 1,2...N) ΠΗ, £4 
Æ CNLSE(3.34) 有 三 孤子 解 ， 
BO em + BD em + BY ens + entnitmt 4 emnes 十 em 起 十 1 二 681 
= ο DNI LEM MI ELE DRIED a 
qeu ΕΠ + eistná-m-6s; 十 emtns+mte 十 ραπ] Πα 101} 4 eren -m 06; 


D 
emtngtm+e; 十 emtnitmtngtnsteis 十 em+nitm-+ngtns+ee; 十 rnm "mega; 
----------ᾱ-----..---ᾱ---τ- τ τν ων νο ος 


D 
(j =1,2,--- ,.N) (3.35) 


其 中 ， 


D 一 1 十 οι η ΓΒ + gn Ra + enitn Bs + entn+ro 
pentiti 4 eMtns+r20 十 mtn +o 十 emtni+vao 
pes TÉ +o 4 ρῃι Γη] Hmm) Ra 十 omtnjtnst+nj+Rs 


renim Re 十 emtnitmt3+wio 十 menfem γω 
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perti m-kn$ hao 十 orn ense tudo 十 erre Γη ἜΣ eso 


«619 Γη eni twig 十 enimse RS (3.36) 
这 里 ， 
Ti 一 k;(T + ik;£), J ταὶ 1,2,3, e j= "YE Me on 十 a 1 
PRECELSI CV ~ PY Kn) ως (2 UL - pP ma) 
E LIO i ΜΠ 
u (Ra (OP ka — BG 2) ος. (Ra = Es) (BY κα — a) 
(kz + Pas Tk) C (kg + 2^ th) C 
PRECELS CL Tu Pro) NR o — BY κοι) 
(Ko + "E 本 2 ᾿ (Κα + 3 ki + ks) j 
_ (ki — ka) (BP καρ — BY) aa) oft = (kı — ks) (8 ios — BY κιο) 
(ki κο) +k) ᾿ (ki +k3)(k3 +k) C 


(Κα — ki) (ks ~ ki) (ks — ka)(k5 — ki) 
(Ei + kj) (Ka + k5) (Kj ko) (Kt + ka) (K2 + ks) (R5 + kı) 
x [68 (Ray Kaa Kazan +BY (κιρκαι--καρκιι)Η-β0 (κα καρ--κιοκοι)]; 
(Κα — k)(ks — Κι)(Κο — ko)(k3 — ki) 
ταις ki) (ks + Eg) (t + a) (t + ka) (3 + Ka) (Rs + kı) 
x (ef (καικαα--καφκαι)Η-β6) (FasKa1 — 3331) 4-03 ϱ x1 κα--κιακοι)], 
(Κα — ki)(ks — #1) (ks — ka)(k3 — k3) 
ee ο 
o 
( 


eh = 了 -一 


KagKs3—Kogks2) +89 (Kran —h33512)-BU (Ki2h23—K22K13)], 
~ (hr + Fi) (ha + ki) (Fs + ἘΠΕ + KP [ka + k3P| ks + kil? 


x [((1ι K22K33 —K11K23K32)--(K12&23 K31 — 12 KagKa1)--(K132 K32— Καὶ 02 K13)]; 


Km K 
dix us m m= 1,2,3 μπε 
ΓΕΒ’ ko + kg’ 


πι. (ha — ki)(k3 ~ Kr) iras — Kyaka) 
(ka + ki) (er + ka) s + R3) eo + k3)’ 
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πε. (ks — Κι) (8 — kr )(&asfii 一 ΚΙΦΚΑΙ) 
^ (1 + ks) (A + KS) (Ks + kg)’ 

Re 一 (ks — Κο)(Κᾷ — k3)(Ka3K22 — Kaska) 
(Ko + k5)(K5 + ks) (Ka + k3 )(ks + k3) 
wo — (ka — Κι)(Κᾷ — Ki)(nitzs — kiska) 
"(hy + RE) (keg + ka)(ki + 3) (ke + 3)’ 

— Lh = ky — ED) (kakaa Kanne) 
(ki + ο) (9 + ka) (s + k3) (Ro + Κᾷ)᾽ 

= (ka — Ki) (ES — ki) (sata 一 askaz) 
(ka + κ.) (ο + ks) (I + k$) (Ka + k3)’ 


3 


3 
a Y fg" , 
κα = i (i, = 1,2,3). k W APG = 1,283, j-12-,N) X 
Δ(Ν + 1) 个 任意 的 复 常数 。 


(a) b) 


lie 1 15 h | 
vi | 由 
p M M Er? ἡ Λ | n | T 
x RR DUE TEmo 
NGN s NN m mt De s 


38: ΞΑ ΚΗΘ ΔΕΑ EIS AEA (3.29) 的 三 孤立 波 的 强度 (a)|pi|? 演化 图 和 相 
应 的 微 扰 方 程 组 (3.28) 的 解 的 模 方 (Ὁ) E? 随时 间 的 演化 图 . 参数 Κι = 1.2 + 0.51, 
ko = 1.2 — 0.51, kg = 1.2— i, € = 0.02, cl = cp = an = an = ax = A = 1, 
By) = By = 8P = By” = BY? = By = Bf? = BY) = BP =1. 
Fig.3-8: Evolution of intensity profile (α)/ρι|2 of three-soliton solutions of CNLSE for 
N=3 and the evolution of intensity profile (b)|E;|? of solution of perturbed CNLSE 
with € = 0.02. The parameters are Κι = 1.2+0.5%, kz = 1.2—0.5i, kg = 1.2 — i, 


= C2 = O11 = 03? = 051 = AQ? = 1, 
go ER gc E a? = a) M p A gp 2 a = ge - a?) =1. 
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3-8(a) SH T ZDERE FRERES I BAN — SUE RIS RE UU 
图 (由 于 [ρ 和 | 到 的 演化 情形 与 |. 和 E.P? SERA, AE 
述 )。 相 应 的 微 扰 CNLSE(3.28) 的 解 的 模 方 的 演化 如 图 3-8(b)。 取 解 (3.33) 中 
的 p, WA (3.35)， 一 阶 修正 q Ha = pr 给 定 。 在 微 扰 项 —ieE 的 作用 下 
三 个 孤子 的 振幅 和 人 能量 都 将 随 着 x 的 增加 而 逐渐 衰减 ， 而 且 微 扰 参量 e RX, 
衰减 得 越 快 。 


(=) (2+1) 维 非 线性 薛 定 刘 方 程 
(2+1) EFRR ESITE 


| 2 o æ 
iu, + Au + |u|*u = 0, A= (gr + ay) (3.37) 


表示 激光 束 在 含有 Kerr FRE AP EMH RARR, Rhu 为 电场 包 
5 ο 为 波 传播 方向 的 轴 向 距离 ，Al 为 横向 拉 普 拉 斯 算 符 ， 它 的 作用 是 产生 
横向 波 矢 kL， 导致 术 射 效应 ，(z,y) 表示 横向 平面 的 坐标 。 等 号 左 端 第 三 项 是 
非 线性 项 。 我 们 知道 ， 该 模型 不 能 用 IST 精确 求解 。 且 (2+1) 维 NLSE 的 解 
可 能 呈现 为 一 有 限时 间 爆 发 或 一 “ 自 聚 焦 ” 的 奇异 性 。 通 常 ， 此 奇异 性 在 对 称 
坐标 下 研究 ， 


Ou 10 δι 2... - = 2 2 
ix alg uud =0, u-wnz, r= y (z? +°). (3.38) 


3,78 9 X (Zakharov) 等 人 发 现 ， 当 选取 合适 的 初始 条 件 w(7,0)， 在 能 量 小 于 零 
时 会 导致 惯量 矩 在 有 限时 间 内 消失 (在 守恒 律 意义 下 )， 结 果 形成 声 缩 (爆发 )， 
这 一 方面 的 详细 内 容 可 以 参见 文献 [91] 及 相关 的 参考 文献 。 

在 研究 二 维 NLSE 时 ， 大 多 数 的 处 理 方法 是 作 变 量 代 换 将 它们 化 为 一 维 情 
况 ， 这 在 数学 上 是 可 以 的 ， 但 在 物理 学 上 有 时 是 不 允许 的 。W.H.Hui,，M .Tajiri, 
A.Nakamura 等 人 却 采 取 了 一 种 简单 解法 对 (2+1) 维 NLSE 进行 了 研究 并 得 到 
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了 它 的 如 下 形式 的 准 孤 立 子 解 92; 
V2wsech za -0t — (9), 
ue strep Sig 20, τα 20) (3.39) 


V5tanh FHE- at ~ à). 
(Zi <0, qa <9) 
其 中 w, ο 为 任意 常数 。 根 据 文献 [91]， 我 们 知道 方程 (3.37) 有 径 向 对 称 的 波 
导 解 : 
u(z,r) = εὖ R(r), 
其 中 OR 满足 非 线性 边界 值 问题 : 


A, R-R+R=0, R(0)=0, Δι el 4 


ro 
该 常 微分 方程 有 很 多 组 解 ， 其 中 最 吸引 人 的 是 Townes 孤子 解 ， 具 体 可 参见 文 
献 93]。 f 

最 近 ， 由 于 Bolti 等 人 的 开拓 性 工作 ，(2+1) HE NLSE 的 dromion 解 也 引 
起 了 人 们 的 关注 : 阮 航 宇 等 人 采用 分 离 变 量 方法 研究 了 (2+1) 维 非 NLSE 的 局 
域 结构 。 他 们 认为 只 要 选择 合适 的 任意 函数 ， 局 域 解 可 以 是 dromion 解 ， 环 孤 
子 解 ， 呼 吸 子 和 瞬 子 等 94, Ad, MOAT ALA Painlevé 性 质 ， 得 
到 了 它 的 无 穷 多 个 对 称 及 它 的 精确 解 (99, 


1 近似 解 
这 一 节 中 ， 我 们 将 考虑 含有 时 间 色 散 的 (2+1) 维 NLSE 
iUz + Uga + Uyy + Juf u + eun = 0. (3.40) 


它 是 描述 在 含有 瞬时 Ker 非 线 性 介质 中 传输 的 超 短 激光 脉冲 的 基本 模 
型 96, 91。 其 中 
oj 
Το : 
AP, το 表示 初始 脉冲 的 半径 ， 波 数 加 = ωπο(ω)/ο, no 为 线性 折射 率 ，c 
和 了 分 别 代表 光速 和 脉冲 持续 时 间 。 若 小 参量 。> 0， 为 正常 色散 ， 若 e < 0， 


则 为 反常 色散 。 当 e = 0， 方程 (3.40) 退化 为 (24-1) 维 NLSE(3.37)。 
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ARORA, HR (3.40) 中 u 展开 成 : 
u = e Du (E, nr, €) 


= 9 (uy(£ n, C,T) + eu(é, m6,7) + O(®)), (3.41) 
其 中 
a=a(z,y,z,t), b= b(z,y,z,t), € = E(z,y,z,t,€), 
n=n(z,y,2,t,€), Ç =C(z,y,2,t,€), T = (y, tie) (3.42) 
均 为 待定 的 实 函 数 。 并 且 £m Cr 有 以 下 的 性 质 : 


{ένη, Cr} S (s, zt], (3.43) 


使 得 e 一 OW, uo 为 无 微 扰 的 (2+1) Æ NLSE 的 准确 解 。 
将 式 (3.41) 代入 微 扰 方程 (3.40)， 得 到 


2ea,, 2,, » 


fuor, + Uogéeés? + osi! + e" ug, uo 
-FieuxcC, + εὐιεεξα᾽ + euin, + e" (euo n" + 2eugu1uo*) 
+e{(i(a, + ib) + (Ger + ἴδια) + (ανν + ibyy)) U0 
«is + 2(a, + ib,)6, + See iyu 


η Nex +N 
(77 + 2(a, + iby) my + ———™ tron 


,Ty 


+(i= + Tes + Tm) 


1 
< F Uor + "s + Cyy toc 


ον ο ος 
ἜΖυρτητψην + 2ρηςΏν(ν) 一 UorrTe }= O(e). (3.44) 
由 于 £n, Gr 有 如 (3.43) 式 的 性 质 ， 我 们 认为 ἐν, ξεν ξεν ἔσαν Πεν Ms Nes Nyys (αν 
Cys Cty τε, Ty，Tz 均 为 e 的 一 阶 项 。 令 近似 方程 (3.44) 中 e 的 相同 次 寡 的 系数 为 
F, 得 


ituocGe + Uogeéz + tomy + €" ug uo" = 0. 
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ieur C, 十 euatetEa2 + CUinn ty? + e" (euo αι” + 2ecuotituo*) 
十 ef (i(az + ibz) + (azz + ibre) + (Qyy + ibyy))uo 
ας 十 2(az + ib,)&, + ἔος + bw + w )uoe 


Ne 


Hie + 2(αν + iby)ny + Nez + Ty 


€ 


Juon 


Tz 


HOZ + yu, + (Can + Gn) 
1 
p (2ύοε-ξατα + 2υηεηξαηα + Ζύοεεξα(α + 2ιοςηξγην 
+2uornTyNy + Ζυρηηγέν) 一 Worr } = O(€ ). (3.45) 
GR, FE uo ETER HRA, PUREE, η, ς 应 满足 以 下 关系 式 ; 
G = ge & =€”, Ty = e^, (3.46) 
而 且 复 函数 uy 为 以 下 无 微 扰 的 (2+1) 4 NLSE 
foc + ogg + Uon + uol uo = 0. (3.47) 


的 精确 解 。 另 外 根据 展开 式 (3.41) 49, αι AAE (8.47) 的 对 称 ， 即 是 它 的 线 
性 化 方程 
fusc + Urge + Uin 十 2upugt + wut = 0 (3.48) 
的 解 。 它 的 最 简单 的 三 个 解 的 形式 可 以 表示 成 1951, 
Ui = Ugg, Ui = Uy, Ul = Um. (3.49) 
给 定 一 个 非 平凡 解 uo, ΠΊΕ (3.45) 和 (3.48) 在 任意 空间 z, y, z 和 时 间 
t 都 是 相 容 的 ， 故 方程 (3.45) 中 的 其 他 项 都 应 等 于 零 ， 即 


i(az + ibz) + (Ger + ἱδσα) + (ayy + iby) = 0, 


i& + 2(a, + ibs )& + Sea + bs t Sw -0 


= 0, 


这 + 2(a, + iby)ny 十 “ee τ y t Duy 


=) (241) 维 非 线性 薛 定 诅 方 程 51 ， 


ÍT, + Tex + Tyy = 0, ἔχη: + ἔγην = 0, 


τὲ =0, Con + Cyy = 0, Exe = 0, 


£T, =0, n6, =0, ηντν = 0. (3.50) 
487538 (3.46) 和 方程 组 (3.50) 得 
a=hz, b= -grm, 


£ = ze, n= ye, 
erezh -1 


2eh 
其 中 hE hlt) 为 关于 时 间 + 的 任意 函数 ，C 为 任意 常数 。 这 样 微 扰 方程 (3.40) 
的 近似 解 便 可 以 写成 : 


T=C, (= 


(3.51) 


hy 224: 2 
u = eO ED e, y, z, t) nlia y zt), ζ(α, y, zt), 7(2,y, 2,t)) 


-Feuy(£(z, y, 2, t), (x, y, 2, ϐ), C(x, y, 2,t), T(z, y, 2,¢))). (3.52) 


2 微 扰 对 孤子 的 影响 


与 (1+1) 维 NLSE 情况 所 不 同 的 是 ， 含 微 扰 项 -eun 的 (2+1) Æ NLSE 的 
微 扰 解 (3.52) 的 振幅 的 变化 指数 不 仅 与 z 有 关 ， 还 与 一 个 关于 时 间 t 的 任意 函 
数 hh 有关。 选取 不 同 的 h 形式 ， 近 似 解 的 形式 会 发 生 较 大 的 改变 。 而 且 ， 不同 
的 h 形式 也 会 对 近似 解 中 的 其 他 参量 产生 影响 。 所 以 ， 在 (2+1) 维 情况 下 ， 微 
扰 项 对 无 微 扰 方 程 的 解 所 产生 的 影响 要 比 (1+1) 维 的 情形 要 复杂 得 多 。 

简单 地 ， 如 取 h(t) 为 常数 ， 则 近似 解 u 的 振幅 只 与 z 有 关 ， 并 随 着 z 的 
变化 而 指数 增加 (z > 0) RER (z < 0)。 若 取 A(t) = t, 根据 (3.39) 式 中 的 情 
形 1, 方程 (3.40) HERA: 


w ,2exte’* — 2ceyte*t 


uo(z, y, z,t) = 2V/ Qwete~**** gech| Tow geet 7 — 6o] 


j er? e2ezt ey2te2est —. 
xd μμ“... (3.59) 
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经 过 简单 的 计算 会 发 现 ， 当 e 一 0， 解 (3.53) 将 退化 回 (2--1) 维 NLSE 的 准 孤 
子 解 。 由 于 它 的 一 阶 修正 更 复杂 ， 这 里 我 们 就 不 再 详细 写 出 。 


(四 ) 小 结 


孤子 及 其 微 扰 问题 的 研究 是 近年 来 非 线 性 领域 一 个 非常 活跃 ， 并 且 硕 果 黑 
累 的 一 个 分 支 。 但 是 一 直 以 来 ， 在 处 理 高 维 的 和 看 合 非 线性 演化 方程 的 微 扰 问 
题 时 ， 现 有 的 微 扰 理论 都 遇 到 了 很 大 的 困难 。 另 一 方面 ， 由 于 高 维 的 和 耦合 非 
线性 演化 方程 已 经 成 功 地 描述 了 实际 问题 中 的 很 多 非 线 性 现象 ， 所 以 又 一 直 受 
到 研究 者 的 广泛 关注 。 

这 一 章 中 ， 我 们 将 直接 微 扰 方 法 成 功 地 应 用 于 含有 线性 损耗 项 的 耦合 非 线 
HEE AEA (2+1) 维 非 线 性 薛 定 廖 方程。 通过 这 种 直接 微 扰 方法 ， 我 们 可 
以 得 到 包 斤 零 阶 近似 和 一 阶 修正 的 可 靠 的 近似 解 。 在 此 基础 上 ， 我 们 还 就 微 扰 
对 孤子 的 各 个 参数 产生 的 影响 进行 了 简单 分 析 : 在 线性 损耗 的 作用 下 ， 孤 子 的 
振幅 随 着 时 间 的 增加 逐渐 衰减 直至 0， 而 孤子 的 传播 速度 并 不 发 生 改 变 。 特 别 
Hh, FAS SESE ES ED MILLER NERS AA MRE RES 
RR. DA- PEKEREN ATAA, BTCA RARE 
定 证 方程 进行 直接 模拟 。 结 果 显示 ， 在 e 值 的 相当 大 的 范围 内 ， 数 值 解 与 解析 
解 吻合 得 很 好 。 但 是 ， 从 解析 解 与 数值 解 的 比较 中 我 们 也 发 现 了 一 个 问题 : EC 
子 的 碰撞 时 刻 并 不 相同 ， 而 且 碰 撞 时 刻 的 孤子 的 振幅 符合 得 也 并 不 是 很 理想 。 
这 是 数值 模拟 过 程 不 够 精确 还 是 直接 微 扰 方 法 在 处 理 多 孤子 碰撞 过 程 时 出 现 了 
问题 ， 还 有 待 于 进一步 研究 。 
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对 含有 修正 项 的 非 线性 演化 方程 的 求解 ， 目 前 的 主要 手段 还 是 停留 在 数值 
方法 。 然 而 当 修 正 项 可 以 看 作 小 量 时 ， 发 展 有 效 的 微 扰 方法 就 成 为 必要 ， 而 且 
各 种 不 同 的 微 扰 方 法 也 是 解决 实际 问题 所 不 可 缺少 的 。 非 线性 芒 定 户 方 程 作为 
现代 科学 中 最 重要 的 也 是 最 普遍 的 非 线性 模型 之 一 ， 关 于 它 的 微 扰 问题 的 研究 
也 越 来 越 受到 物理 学 家 和 数学 家 们 的 广泛 关注 。 到 目前 为 止 ， 人 们 也 已 经 提出 
了 很 多 微 扰 方法 来 研究 非 线性 薛 定 户 方 程 的 微 扰 问 题 ， 包 括 逆 散 射 徽 扰 方 法 ， 
直接 微 扰 方法 ， 修 正 守恒 律 微 扰 方 法 ， 奇 点 微 扰 方法 ， 线 性 微 扰 方 法 ， 拉 格 朗 
日 微 拢 方法 等 等 。 最 近 ， 楼 在 前 人 的 基础 上 发 展 了 非 线性 薛 定 刘 方 程 的 一 种 新 
的 直接 微 扰 方 法 。 本 学 位 论文 中 ， 我 们 系统 地 研究 了 楼 提出 的 新 的 直接 微 扰 方 
法 ， 并 将 其 应 用 于 微 扰 耦 合 非 线 性 薛 定 诅 方程 和 (2+1) SEARS HERE EIS FE. 

较 之 其 他 的 微 扰 方法 楼 的 直接 微 扰 方 法 具有 一 定 的 优越 性 ，(1) 完全 摆脱 
了 对 逆 散 射 方法 的 依赖 。 也 就 是 在 这 种 直接 微 扰 方 法 中 ， 不 要 求 未 受 微 扰 影 
响 的 非 线 性 薛 定 廖 方程 能 够 用 逆 散 射 法 精确 求解 ，(2) 思路 直接 ， 容 易 理 解 和 
接受 ，(3) 同时 得 到 各 种 不 同类 型 的 微 扰 解 ， 包 括 单 孤子 ， 多 孤子 ， 周 期 解 等 
等 ， 而 不 需要 采用 不 同 的 微 扰 方 法 来 分 别处 理 它们 的 微 扰 问题 ， 尤 其 象 一 直 受 
到 人 们 关注 的 多 孤子 微 扰 问题 ， 楼 的 直接 微 扰 方 法 也 能 很 好 地 处 理 。(4) 用 这 
种 直接 微 扰 方法 得 到 的 近似 解 可 以 以 无 微 扰 的 非 线 性 苹 定 户 方 程 的 任意 一 个 精 
确 解 作为 它 的 零 阶 解 ， 而 且 其 一 阶 修正 可 以 由 它们 的 无 穷 多 个 对 称 中 的 任意 一 
个 给 定 。 所 以 ， 运 用 这 种 直接 微 扰 方 法 得 到 的 近似 解 比 用 通常 的 微 扰 方法 得 到 
的 近似 解 要 丰富 得 多 。 通 过 将 直接 微 扰 方法 运用 到 一 个 可 积 的 含 修正 项 的 非 线 
性 薛 定 刘 方 程 ， 我 们 对 它 进 行 了 可 车 性 分 析 。 从 近似 解 与 精确 解 的 比较 中 发 
现 ， 在 相当 大 的 e 值 范围 内 ， 近 似 解 与 精确 解 符合 得 很 好 。 于 是 在 论文 的 第 三 
章 中 ， 我 们 将 直接 微 扰 方法 应 用 到 微 扰 非 线性 薛 定 廖 方程 ， 并 得 到 了 它们 的 可 
人 靠 的 近似 解 。 楼 的 直接 微 扰 方法 还 能 够 给 出 孤子 参数 的 解析 表达 式 ， 且 与 数值 
结果 符合 得 非常 好 。 在 线性 损耗 的 作用 下 ， 孤 子 的 振幅 随 着 时 间 的 增加 逐渐 误 
减 直至 0， 而 孤子 的 传播 速度 并 不 发 生 改 变 。 特 别 地 ， 耦 合 非 线 性 薛 定 户 方程 
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的 双 孤 立波 的 能 量 不 会 因为 微 扰 的 作用 而 发 生 转 移 或 交换 。 

诚然 ， 本 文 所 做 的 工作 仅仅 只 是 一 个 开始 ， 还 有 许多 的 问题 有 待 进一步 的 
发 展 和 完善 : 

1. 虽然 我 们 成 功 地 处 理 了 几 个 具体 的 含 不 同 微 扰 的 非 线性 薛 定 刘 方程 的 例 
子 ， 但 是 目前 我 们 还 无 法 将 直接 微 扰 方 法 运用 到 含 其 他 形式 的 微 扰 的 非 线 性 薛 
BA, W icus u, ielultu, ieju[^u, 等 等 ， 其 中 的 难点 在 于 我 们 无 法 找到 
满足 相应 的 相 容 性 条 件 的 解 。 

2. 至 于 如 何 将 该 直接 微 扰 方 法 推广 到 其 他 的 微 扰 非 线性 演化 方程 ， 如 微 
ti KdV 方程 、sine-Gordon 方程 、Burgers 方程 ， 甚 至 是 微 扰 的 离散 非 线性 演 
化 方程 等 等 ， 这 些 问 题 都 有 待 于 进一步 讨论 。 

3. 从 第 三 章 处 理 看 合 非 线性 苹 定 计 方 程 的 微 扰 问 题 的 过 程 中 我 们 发 现 ， 在 
双 孤 子 的 碰撞 时 刻 ， 利 用 直接 微 扰 方 法 得 到 的 解析 解 与 数值 结果 出 现 了 明显 的 
差别 ， 是 在 现 有 的 框架 内 进一步 完善 还 是 寻求 新 的 方法 来 解决 ， 尚 待 进一步 深 
入 研究 。 
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分 步 傅立叶 方法 


分 步 傅立叶 方法 (Split-step Fourier Method) 最 早 是 在 1973 年 开始 应 用 
的 98， 这 种 方法 已 经 被 证 明 是 相同 精度 下 数值 求解 非 线性 花 定 谓 方 程 最 快 
的 方法 ， 部 分 原因 是 它 采 用 了 快速 傅立叶 变换 算法 (Fast Fourier Transform 
Algorithm)。 由 于 它 比 大 多 数 有 限 差分 法 见效 快 ， 已 被 广泛 地 用 来 处 理 非 线性 
色散 介质 地 脉冲 传输 问题 的 非 线性 薛 定 请 方程 89, 100。 下 面 我 们 就 分 步 傅 立 
叶 方 法 作 一 简单 的 介绍 。 

非 线性 薛 定 谓 方程 的 基本 形式 为 : 


iu, = ug, η (A.1) 


其 中 必 是 复 值 的 未 知 函 数 。 
设 当 |z| oo t, juj +0, MARR EEA (A.1) 的 初 值 问题 有 单 
孤立 子 的 精确 解 


u(a, t) = 2ge ρε: (ο) mt dol sech (Ina — 8ENt 一 zo) (A.2) 


其 中 En poro HARK. BR, (Al) 的 初 值 也 可 以 由 两 个 孤立 子 甚至 多 个 
孤立 子 的 精确 解 给 定 。 

为 了 叙述 的 简洁 ， 我 们 考虑 NLS 方程 (A.1) 的 周期 初 值 问题 。 初 始 条 件 和 
边界 条 件 是 


u(x, t)li=0 = uo(z) 
u(x + 2π, 8) = u(z,t). 


将 区 间 [0, 27] 分 为 J 等 份 ， Zo = 0, zı = Az μα; = 27, Ar= x, 
这 种 方法 对 方程 (A.1) 分 步 求解 
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(1) 首先 考虑 非 线性 部 分 
iuf) = Qu Put, (A3) 
ul) [150 = Ulta = Uo(Z). 
它 有 精确 的 单 波 解 
u(x,t) = ug(z) - e 2iluo(z)|?-t (A.4) 
(2) 再 考虑 线性 部 分 
iuf = ull), (Α.5) 
用 Fourier 变换 求解 ， 在 Fourier 空间 中 有 解 
al = a" (gyg^t. (A.6) 


(3) 将 (A.4) 的 解 u(x, t) 变换 到 Fourier 空间 为 如) (以 它 作为 (A.5) 
中 的 初 值 
a (t+ At) = (te δ', 
(4) 将 GU (t e. At) 变换 回 x 的 空间 中 得 到 近似 解 
aj(t+ At) = F (00D (t+ ΔΘ} = Fe ^t. F(uD(z,t)), — (AT) 


这 个 解 能 够 用 FFT 标准 程序 ， 很 快 求 得 。 
这 种 方法 对 时 间 t 是 二 阶 截断 误差 。 而 对 空间 α 的 截断 误差 为 (Az), 
其 中 M 取决 于 解 的 光滑 性 。 当 解 无 穷 次 可 微 时 ，M 可 以 取 任何 正 值 。 
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